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In the following paper several problems are considered, wliich were alluded to in my 
two previous papers on this subject.''^' 

The paper is divided into sections which deal with the problems referred to in the 
table of contents. It was foimd advantageous to throw the several investigations 
together, because their separation would have entailed a good deal of repetition, 
and one system of notation now serves throughout. 

It has, of course, been impossible to render the mathematical parts entirely inde- 
pendent of the previous papers, to which I shall accordingly have occasion to make a 
good many references. 

As the whole inquiry is directed by considerations of apphcability to the earth, I 
shall retain the convenient phraseology afforded by speaking of the tidally distorted 
spheroid as the earth, and of the disturbing body as the moon. 

It is probable that but few readers will care to go through the somewhat complex 
arguments and analysis by which the conclusions are supported, and therefore in the 
fourth part a summary of results is given, together with some discussion of their 
physical applicability to the case of the earth, 

I. Secular distortion of the spheroid^ and certain tides of the second order, 

* 

In considering the tides of a viscous spheroid, it was supposed that the tidal protu- 
berances might be considered as the excess and deficiency of matter above and below 

* ** On the Bodily Tides of Viscous and Semi-elastic SpheroidSj and on the Ocean Tides upon a Yielding 
Kucleus," Phil. Trans., 1879, Part I., and — 

*' On the Precession of a Viscous Spheroid, and on the remote History of the Earth," immediately pre- 
ceding the present paper. They will be referred to hereafter as " Tides " and '^ Precession " respectively." 
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the mean sphere — or more strictly the mean spheroid of revokition which represents 
the average shape of the earth. The spheroid was endued with the power of gravi- 
tation, and it was shown that the action of the spheroid on its own tides might be 
found approximately by considering the state of flow in the mean sphere caused by 
the attraction of the protuberances, and also by supposing the action of the protu- 
berances on the sphere to be normal thereto, and to consist, in fact, merely of the 
weight (either positive or negative) of the protuberances. 

Thus if a be the mean radius of the sphere, w its density, g mean gravity at the 
surface, and r=:a'^ai the equation to the tidal protuberance, where cri is a surface 
harmonic of order i, the potential per unit volume of the protuberance in the interior 



r\ ^ 



of the sphere is -^r-^\ -) cri, and the sphere is subjected to a normal traction per unit 

area of surface equal to -—gwcri. 

It was also shown that these two actions might be compounded by considering the 
interior of the sphere (now free of gravitation) to be under the action of a potential 

2(i-l) /rV 

This expression therefore gave the effective potential when the sphere was treated 
as devoid of gravitational power. 

It was remarked''' that, strictly speaking, there is tangential action between the pro- 
tuberance and the surface of the sphere. And latert it was stated that the action of 
an external tide-generating body on the lagging tides was not such as to form a 
rigorously equilibrating system of forces. The effects of this non-equilibration, in as 
far as it modifies the rotation of the spheroid as a whole, were considered in the paper 
on '^ Precession. ^^ 

It is easy to see from general considerations that these previously neglected tangential 
stresses on the surface of the sphere, together with the effects of inertia due to the 
secular retardation of the earth's rotation (produced by the non-equilibrating forces), 
must cause a secular distortion of the spheroid. 

This distortion T now propose to investigate. 

In order to avoid unnecessary complication, the tides will be supposed to be raised 
by a single disturbing body or moon moving in the plane of the earth's equator. 

Let Q^=za-\-o- be the equation to the bounding surface of the tidally-distorted earth, 
where cr is a surface harmonic of the second order. 

I shall now consider how the equilibrium is maintained of the layer of matter cr, as 
acted on by the attraction of the spheroid and under the influence of an external dis- 
turbing potential V, which is a solid harmonic of the second degree of the coordinates 
of points within the sphere. J The object to be attained is the evaluation of the stresses 

*' ''Tides/' Section 2. 
t " Tides," Section 5. 
I A parallel invesbigation would be applicable, where o- and V are of any orders. 
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tangential to the surface of the sphere, which are exercised by the layer cr on the 
sphere. 

Let 9, (f) be the colatitude and longitude of a point in the layer. Then consider a 
prismatic element bounded by the two cones $, $-{-89, and by the two planes cj), ^+8^. 

The radial faces of this prism are acted on by the pressures and tangential stresses 
communicated by the four contiguous prisms. Bat the tangential stresses on these 
faces only arise from the fact that contiguous prisms are solicited by slightly different 
forces, and therefore the action of the four prisms, surrounding the prism in question, 
must be pxincipally pressure. I therefore propose to consider that the prism resists 
the tendency of the impressed forces to move tangentially along the surface of the 
sphere, by means of hydrostatic pressures on its four radial faces, and by a tangential 
stress across its base. 

This approximation by which the whole of the tangential stress is thrown on to the 
base, is clearly such as slightly to accentuate, as it were, the distribution of the 
tangential stresses on the surface of the sphere, by which the equilibrium of the layer 
cr is maintained. For consider the following special case : — Suppose cr to be a surface of 
revolution, and V to be such that only a single small circle of latitude is solicited by a 
tangential force everywhere perpendicular to the meridian. Then it is obvious that, 
strictly speaking, the elements lying a short way north and south of the small circle 
would tend to be carried with it, and the tangential stress on the sphere would be a 
maximum along the small circle, and would gradually die away to the north and 
south. In the approximate method, however, which it is proposed to use, such an 
application of external force would be deemed to cause no tangential stress to the 
surface of the sphere to the north and south of the small circle acted on. This special 
case is clearly a great exaggeration of what holds in our problem, because it postulates 
a finite difference of disturbing force between elements infinitely near to one another. 

We will first find what are the hydrostatic pressures transmitted by the four 
prisms contiguous to the one we are considering. 

Let p be the hydrostatic pressure at the point r, 9, (\> of the layer cr. Then if we 
neglect the variations of gravity due to the layer cr and to V, p is entirely due to the 
attraction of the mean sphere of radius a. 

The mean pressure on the radial faces at the point in question is ^gwcr ; where cr 
is negative the pressures are of course tractions. 

We will first resolve along the meridian. 

The excess of the pressure acting on the face 9'\-h9 over that on the face 9 (whose 
area is era sin 9h^) is 

d . d , 

-T^[^g't(;cr. era sin 9Scj>]89y or ^gi^ctj2){^'^ ^i^ ^)8^S</), 

and it acts towards the pole. 

The resolved part of the pressures on the faces </)+S^ and </> (whose area is cra89) 
along the meridian is 
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{^gwcr) (craSO) (cos OS(j)) or ^gwacr' cos OSOScj), 

and it acts towards the equator. 

Hence the whole force due to pressure on the element resolved along the meridian 
towards the equator is 

^gwaSOScl){cr^ cos ^~";^(o^^ sin 6)), or -—gwaSdScj) sin Oaj^. 

But the mass of the elementary prism Sm=iva^ sin OS0S(f),(T, 

Hence the meridional force due to pressure is — Smy^, 

We will next resolve the pressures perpendicular to the meridian. 
The excess of pressure on the face </)+S</> over that on the face </> (whose area is 
craBO), measured in the direction of </> increasing, is 

Hence the force due to pressure perpendicular to the meridian is — -Sm-r— 7 tT' 

We have now to consider the impressed forces on the element. 

Since cr is a surface harmonic of the second degree, the potential of the layer of 

matter cr at an external point is fgcrl-] . Therefore the forces along and perpen- 
dicular to the meridian on a particle of mass 8m, just outside the layer a but infinitely 

near the prismatic element, are i-Sm-r;: and i-Sm -r—^ -7-5 a^nd these are also the forces 
^ ^a dO ^a ^mddcj) 

acting on the element 8m due to the attraction of the rest of the layer cr. 

Lastly, the forces due to the external potential V are clearly Sm--j^ and 



dd 



dm — r-^jr. 
a sm 6 d(p 

Then collecting results we get for the forces due both to pressure and attraction, 
along the meridian towards the equator 



8m 



g dcr ^0 dcr dY 
adO ^adO adO 



^"^^e^^ -''^""^ 



and perpendicular to the meridian, in the direction of ^ increasing, 



hm 



a d(T , 3(7 d<j , dN 



^'^T5::ra^i^ ^^9^)^ 



a sin d d(f> 5a sin d^ a sin 6d(f>_\ a sin dcj) 

Henceforward — will be written g, as in the previous papers. 
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Now these are the forces on the element which must be balanced by the tangential 
stresses across the base of the prismatic element. 

It follows from the above formulas that the tangential stresses communicated by the 
layer cr to the surface of the sphere are those due to a potential V—gacr acting on 
the layer a, 
V 



If 



cr: 



^ct 



there is no tangential stress. But this is the condition that cr should be 



the equilibrium tidal spheroid due to V, so that the result fulfils the condition that if 
o" be the equihbrium tidal spheroid of V there is no tendency to distort the spheroid 
further ; this obviously ought to be the case. 

In the problem before us, however, cr does not fulfil this condition, and therefore 
there is tangential stress across the base of each prismatic element tending to d.istort 
the sphere. 

Suppose V=r^S where S is a surface harmonic. 

Then at the surface V=:a^S. If hm be the mass of a |)rism cut out of the layer cr, 
which stands on unit area as base, then hm=^iD(T, 

Therefore the tangential stresses per unit area communicated to the sphere are 



d 



iva^~ jz {B — Q-) along the meridian 



a dO 



and 



> 



(Tx 



%va^- - — - -77 (S — Q:") perpendicular to the meridian 
a^md d(p ^ ^cr ^ ^ 



(1) 



Besides these tangential stresses there is a small radial stress over and above the 
radial traction— p^M^cr, which was taken into account in forming the tidal theory. But 
we remark that the part of this stress, which is periodic in time, will cause a very 
small tide of the second order, and the part which is non-periodic will cause a very 
small permanent modification of the figure of the sphere. But these effects are so 
minute as not to be worth investigating. 

We will now apply these results to the tidal problem. 
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Let XYZ (fig. 1) be rectangular axes fixed in the eartli, Z being the axis of 
rotation and XZ the plane from which longitudes are measured. 

Let M be the projection of the moon on the equator, and let m be the earth's 
angular velocity of rotation relatively to the moon. 

Let A be the major axis of the tidal ellipsoid. 

Let A'K=o)t^ where t is the time, and let MA=e. 

Let m be the moons mass measured astronomically, and c her distance, and 



3^ 



Then according to the usual formula, the moon's tide-generating potential is 



Tr^[sin^ cos^ {(fy—cot — c)-— ^], 



which may be written 



|Tr^(-|-— cos^ 6)'-^^Tr^ sin^ cos 2(^— coi^ — c). 

The former of these terms is not a function of the time, and its effect is to cause a 
permanent small increase of ellipticity of figure of the earth, which may be neglected. 
We are thus left with 

^TT^ sin^ cos 2(<^ — ft)^ — e) 
as the true tide-generating potential. 

1 uUCO 

Now if tan 2€=— — ; where v is the coefficient of viscosity of the spheroid, then by 
the theory of the paper*on ^^ Tides, ^^ such a potential will raise a tide expressed by 



- =1- cos 2e sin^ cos 2(^ — a)^)"'" 
a g 



^9^ 

» « » • • • \ iJ I 



ihen II we put 



S=-|Tsin^ ^cos 2(^— o)^ — e) . ....... (3) 



and 



S— g~=-|T sin 2e sin^ sin 2{<j)-—o)t) . . . 



d 



cl 



-f S— 'g- ]=T sin 2c sin cos sin 2(^— co^) 



• (4) 



1 d 



mi6 cl(jy 



S— •g-^j=T sin 2e sin cos 2{(f>-—a)t), 



Multiplying these by 'i/ja^-, we find from (1) the tangential stresses communicated 
by the layer <t to the sphere. 



^ " TideSj" Secfcion 5. 



^ 
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They are 



rr-i 



tt'Ci^l—- sin 4e sin^ 6 cos 6 sin 4{(j>—o)t) along the meridian, 
S 



and 



rr-'a 



iva^- sin ie sin^ ^(1+ cos 4((^— wi^)) perpendicnlar to the meridian. 
S 



These stresses of course vanish when e is zero, that is to say when the spheroid 
is perfectly fluid. 

In as far as they involve (ft—mt these expressions are periodic, and the periodic 
parts must correspond with periodic inequalities in the state of flow of the interior of 
the earth. These small tides of the second order have no present interest and 
may be neglected. 

We are left, therefore, with a non-periodic tangential stress per unit area of the 
surface of the sphere perpendicular to the meridian from east to west equal to 



\yja?— sin 4e sin*^ 0, 



The sum of the moments of these stresses about the axis Z constitutes the tidal 
frictional couple J^, which retards the earth's rotation. 
Therefore 



\ivo?- sin 4€ sin^ ^. a sin ^. a^ sin 6d6d(j) 



B 



integrated all over the surface of the sphere, and efiecting the integration we have 



4:7r f. T^ , . 
waP , — sm 4e. 



15 g 

But if C be the earth's moment of inertia, Q-=--^^inva^, 
Therefore 



s_i 



3 



C —f^ sin 4€ . . . (5) 

This expression agrees with that found by a different method in the paper on 
''Precession."''^ 

We may now write the tangential stress on the surface of the sphere as ^iva^~ sin^ 9 ; 

and the components of this stress parallel to the axes X, Y, Z are 



--|ira^77sin^^sin^, +i'^^<^^7rsin^ ^cos ^, . . . . . (6) 

We now have to consider those effects of inertia which equilibrate this system of 
surface forces. 

The couple J^ retards the earth's rotation very nearly as though it were a rigid 

# "Precession," Section 5 (22), wlien {=0, 

4 A 2 



546 



ME. a. H. DARWIN" ON PROBLEMS CONNECTED 



body. Hence the effective force due to inertia on a unit of volume of the interior of 



the earth at a point r, 6, (f) is tor sin -j^^ and it acts in a small circle of latitude from 

west to east. The sum. of the moments of these forces about the axis of Z is of course 
equal to J^, and therefore this bodily force would equihbrate the surface forces found 
in (6), if the earth were rigid. 

The components of the bodily force parallel to the axes are in recta-ngular co- 
ordinates. 

m 



toy- 



c c 



(^) 



The problem, is therefore reduced to that of .finding the state of flow in the interior 
of a viscous sphere, which is subject to a bodily force of which the components are (7) 
and to the surface stresses of which the components are (6). 

Let a, ^, y be the component velocities of flow at the point x, y^ %^ and v the 
coefiicient of viscosity. Then neglecting inertia because the motion is very slow, the 
equations of motion are 



-^^+vV''-a-\-iD^y 
cl/p 



, A-vV^B—iv-zzx 

ay U 







= 






+f V^7 



da , <^^y8 |^% 

dx dy dz 



= Q 



= 



. . . (8) 



_y 



We have to find a solution of these equations^ subject to the condition above stated^ 
as to surface stress. 

Let ol\ ^\ y\ p be functions which satisfy the equations (8) throughout the sphere. 
Then if we put a=ia^a^^ ^=yS^4"A? y^^T'"4"7/? J^=p'+P/? we see that to complete 
the solution we have to find ex^, ^^, y^, j:)^, as determined by the equations 



' dy d^ dx dy 



dz 



dz 



• • (9) 



which they are to satisfy throughout the sphere. They must also satisfy certain 
equations to be found by subtracting from the given surface stresses (6), components 
of surface stress to be calculated from a , ^\ y\ p\ "''' 

We have first to find a\ ^\ y\ p\ 

Conceive the symbols in equations (8) to be accented, and differentiate the first 



* This statement of method is taken from Thomson and Tait's ^ Fat. Pliil.,' § 733. 



WITH THE TIDES OF A VISCOUS SPHEROID. 



547 



three by x, y, z respectively and add them ; then bearing in mind the fourth equation, 
we have v^p'=0, of which ^^'=0 is a sohition. 
Thus the equations to be satisfied become 



V 2a': 



y, V2^=- Y^x, vV=0. 



V C 



V G 



Solutions of these are obviously 






1 !f#,,2, 
V c 



10.. ri''^ y> 



yS'= 



1- 1^ m^% 

V C 



10 .. ri ' •^^ 



7=0 



IV 



10 



V C 



r^ sin 6 sin ^ =tV ~ 77'^'^ ^i^^ ^ cos ^ 



>- 



V C 



. (10) 



These values satisfy the last of (8), viz. : the equation of continuity, and therefore 
together with ^'=0, they form the required values of a^ ^\ y\ p\ 

We have next to compute the surface stresses corresponding to these values. 

Let P, Q, B, S, T, U be the normal and tangential stresses (estimated as is usual 
in the theory of elastic solids) across three planes at right angles at the point x^ y, z. 

Then 



■^ clx \ dz di 



(11) 



Q, R, T, U being found by cyclical changes of symbols. 

Let F, G, H be the component stresses across a plane perpendicular to the radius 
vector T at the point x, y, z ; then 



Pr=P^+Uy+T^ 



> 



(12) 



d 
Substitute in (12) for P, Q, &c., from (U), and put ^'— a^.^r+iS^^Z+y'^.* ^^^^ vj- for 

CO Cv Co rfM 

x-r + y- — \-z~, I hen 

dx ^ dy dz 



Fr 



— px+vi \T~ 1 \a:A — -- 

•^ \ clr / dx 



^, G7"=&c., Hr=&c. 



(13) 



These formulas give the stresses across any of the concentric spherical surfaces. 
In the particular case in hand p'=0, 7^=0, ^=0^ and a! , /3' are homogeneous functions 
of the third degree, hence 
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F= —jiv^^r^ sin 9 sin (b, G=^liv~^-r^ sin 6 cos c^^ 



. 



« * 



(14) 



and at the surface of the sphere r=a. 

Then according to the principles above explained^ we have to find a^, fi^, y^ so that 
they may satisfy 

— -f^ + i;V% =0, &C., &C., 
dx ' 

throughont a sphere, which is subject to surface stresses given by subtracting from (6) 
the surface values of F, G, H in (14), Hence the surface stresses to be satisfied by 
^/' A? 7/' li^ve components 



XJLc 



V) 



"tt \j 



a^'(f — sin^ 6) sin B sin ^, B, 



— 7^«'(l" sin^ B) sin B cos ^, Cg^O, 



These are surfgice harmonics of the third order as they stand, 

Nov/ the solution of Sir W. Thomson's problem of the state of strain of an incom- 
pressible elastic sphere, subject only to surface stress, is applicable to an incompressible 
viscous sphere, mutatis mutandis. His solution"''' shows that a surface stress, of which 
the components are A.;, B^, C^- (surface harmonics of the i^^' order), gives rise to a state 
of flow expressed by 



1 



ct — 



vcc 



i-1 



(a'-r^) cm__^ 1 



_2(2i^-{'l) dx ' i-1 



(^ + 2)r 



2^+1 



d 



_(2'^^^ + l)(2^+l) dx 



{^ 



1 



^i^i^ 



<l>+q 



>'^2i{2i + l) dx ^ ^^ _ 



(15) 



and symmetrical expressions for /3, y. 

Where "i^ and # are auxiliary functions defined by 



%-^ 



t^y 



^(A,rO+|(B.rO+i(C.rO 



^.+1=^ 



dx 

2R3 



dz 



d 

da 



i^p 



,— ?'— 1 



) + ~(B.r-"^) + IXC.^^-^"""^) 



d 



> 



• • (16) 



In our case {=3, and it is easily shown that the auxiliary functions are both zero, 
so that the requ.ired solution is 



/v ■ ■ 1 -2 



r> r, vf — sin^' 0)r^ sin 6 sin 6, B. 



10 



Sv C 



-^(-f — sin^ 0) sin 6 cos (f>, y^ = 0. 



If we add to these the values of a, fi\ y from (10)^ we have as the complete 
solution of the problem, 

* Thomson and Tait's ' Nat. Pliil.,' § 737, 
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a=: — --- "~:^r^ sin^ sin 6, l3-=s- "^r^ sin^ cos 6, y=0 . ... (17) 

8f U ' Of (J 

These values show that the motion is simply cylindrical round the earth's axis, each 
point moving in a small circle of latitude from east to west with a linear velocity 

~ ~r^ sin^ 6, or with an angular velocity about the axis equal to ^ -^r'^ sin^ 6/'' 

In this statement a meridian at the pole is the curve of reference, but it is more 
intelligible to state that each particle moves from west to east with an angular velocity 

about the axis equal to -- ^(a^ — r^ sin^ 6), with reference to a point on the surface at 
the equator. 

The easterly rate of change of the longitude L of any point on the surface in 

colatitude d is therefore - - ~— cos^ d, 

inen smce --r=~ sm 2e cos 2e, and tan 2e=4- t? tnereiore 

'f =it(g cos 2e)'a> cos^ ^ (17') 

This equation gives the rate of change of longitude. The solution is not apph cable 
to the case of perfect fluidity, because the terms introduced by inertia in the equations 
of motion have been neglected ; and if the viscosity be infinitely small, the inertia 
terms are no longer small compared with those introduced by viscosity. 

In order to find the total change of longitude in a given period, it will be more 
convenient to proceed from a different formula. 

Let 7^, /2 be the earth's rotation, and the moon's orbital motion at any time ; and let 

the sufiix to any symbol denote its initial value, also let £=:(--^ 

Then it was shown in the paper on ^^ Precession" that the equation of conservation 
of moment of momentum of the moon-earth system is 



n 



% 



l+/x(l — |^)t (18) 



Where /x is a certain constant, which in the case of the homogeneous earth with the 
present lengths of day and month, is almost exactly equal to 4. 
By differentiation of (18) 

dn d^ /i n\ 

~d-~^'^Ht ^^^> 

'^ The probleni might probably be solved more shortly without using the general solution, but the 
general solution will be required in Part III. 

f "Precession," equation (73), when ^=0 and t' = 0. 
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But the equation of tidal friction is --7=:::—^. Therefore 



dt jjb C% 

Now 

clh loa^ ^ 



at ^v \j 
Therefore 

-TT —— lXt(/r\ ~Z GOd U» « • • . « » . . . { JL t/ ) 

d^ ' ^ Sv ^ ' 

All the quantities on the right-hand side of this equation are constant, and there- 
fore by inteOTation we have for the chano^e of lons^itude 



VXi? 



But since iii^'=-n^'—Sl^^ and tan 2eQ=f • -— ~^, therefore in degrees of arc, 

LJLqU it/ iho."—' aun \ r\ I i* -i \ /J 

__ = /^^Ooa: ^ cot 2€Q(f — 1) COS"^. 

In order to make the numerical results comparable with those in the paper on 
*' Precession/' I will apply this to the particular case which was the subject of the first 
method of integration of that paper. "^ It was there supposed that 6q=:17° 30', and it 
was shown that looking back about 46 million years |^ had fallen from unity to *88. 
Substituting for the various quantities their numerical values, I find that 

~ AL=(f*31 cos^^=19'cos^ B. 

Hence looking back 46 million years, we find the longitude of a point in latitude 
80°, further west by 4f' than at present, and a point in latitude 60°, further west by 
1 4^^ — both being referred to a point on the equator. 

Such a shift is obviously quite insignificant, but in order to see whether this 
screwing motion of the earth's mass could have had any influence on the crushing of 
the surface strata^ it will be well to estimate the amount by which a cubic foot of the 
earth's mass at the surface would have been distorted. 

The motion being referred to the pole, it appears from (17) that a point distant p 

from the axis shifts through r- ^—p^St in the time 8^. There would be no shearing if 

^ ov L' 

# ii Precession/' Section 16. 
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a point distant p+Sp shifted through ~- ~ p^{p-j'Bp)8t ; but this second point does 



shift through - ^{p-^SpfSt 

OV \J 

Hence the amount of shear in unit time is 



1 to 
— ^„ 

Sp 8v C 



i^p-^op) — {p-^op)p 



% 






Therefore at the equator^ at the surface where the shear is greatest^ the shear per 
unit time is 






With the present values of r and cy^ -^f [ - j oj is a shear of Tpj- per annum. 

Hence at the equator a slab one foot thick would have one face displaced with 
reference to the other at the rate of j-q^q cob^ 2€ of an inch in a million years. 

The bearing of these results on the history of the earth will be considered in 
Part IV. 



The next point which will be considered is certain tides of the second order. 

We have hitherto supposed that the tides are superposed upon a sphere ; it is, how- 
ever^ clear that besides the tidal protuberance there is a permanent equatorial pro- 
tuberance. Now this permanent protuberance is by hypothesis not rigidly connected 
with the mean sphere; and, as the attraction of the moon on the equatorial regions 
produces the uniform precession and the fortnightly nutation, it might be (and indeed 
has been) supposed that there would arise a shifting of the surface with reference to 
the interior, and that this change in configuration would cause the earth to rotate 
round a new axis, and so there would follow a geographical shifting of the poles. I 
will now show, however, that the only consequence of the non-rigid attachment of the 
equatorial protuberance to the mean sphere is a series of tides of the second order in 
magnitude, and of higher orders of harmonics than the second. 

For a complete solution of the problem the task before us would be to determine 
what are the additional tangential and normal stresses existing between the protu- 
berant parts and the mean sphere, and then to find the tides and secular distortion (if 
any) to which they give rise. 

The first part of these operations may be done by the same process which has just 
been carried out with reference to the secular distortion due to tidal friction. 

The additional normal stress (in excess of —gwa-, the mean weight of an element of 
the protuberance) can have no part in the precessional and nutational couples, and the 
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remark may be repeated that, that part of it which is non-periodic will only cause a 
minute change in the mean figure of the spheroid which is negligeable, and the part 
which is periodic will cause small tides of about the same magnitude as those caused 
by the tangential stresses. With respect to the tangential stresses, it is d priori 
possible that they may cause a continued distortion of the spheroid, and they will 
cause certain small tides, whose relative importance we have to estimate. 

The expressions for the tangential stresses, which we have found above in (1), are 
not linear, and therefore we must consider the phenomenon in its entirety, and must 
not seek to consider the processional and nutation al effects apart from the tidal effects. 

The whole bodily potential which acts on the earth is that due to the moon (of 
which the full expression is given in equation (3) of '^ Precession^'), together with that 
due to the earth's diurnal rotation (being l-nV^ (^— cos^^)) ; the whole may be called 
r^S. The form of the sru^face cr is that due to the tides and to the non-periodic part 

Ct 11/ 

of the moon s potential, together with that due to rotation — being - ~ (-§-— cos^ 6), 

Now if we form the effective potential a^ (S — g ~ ), which determines the tangential 

stresses between cr and the mean sphere, we shall find that all except periodic terms 
disappear. This is so whether we suppose the earth's axis to be oblique or not to the 
lunar orbit, and also if the sun be supposed to act. 

Then if we differentiate these and form the expressions 

,^ a d I ^ a\ ^^ a d I ^ a 



a d6 \ ^ a/' a sin 0d(p 



a. 



We shall find that there are no non-periodic terms in the expression giving the tan- 
gential stress along the meridian ; and that the only non«periodic terms which exist in 
the expression giving the tangential stress perpendicular to the meridian are precisely 
those whose effects have been already considered as causing secular distortion, and 
which have their maximum effect when the obliquity is zero. 
Hence the whole result must be — 

(1) A very minute change in the permanent or average figure of the globe ; 

(2) The secular distortion already investigated ; 

(3) Small tides of the second order. 

The one question which is of interest is, therefore— Can these small tides be of any 
importance ? 

The sum of the moments of all the tangential stresses which result from the above 
expressions, about a pair of axes in the equator^ one 90° removed from the moon's 
meridian and the other in the moon's meridian, together give rise to the processional 
and nutational couples. 

Hence it follows that part of the tangential stresses form a non-equibrating system of 
forces acting on the sphere's surface. In order to find the distorting effects on the globe, 
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we should, therefore^ liave to eqiiibrate the system by bodily forces arising frora the 
effects of the inertia due to the uniform precession and the fortnightly nutation— just 
as was done above with the tidal friction. This would be an exceedingly laborious 
process ; and although it seems certain that the tides thus raised would be very small, 
yet we are fortunately able to satisfy, ourselves of the fact more rigorously. Certain 
parts of the tangential stresses do form an equibrating system of forces, and these are 
precisely those parts of the stresses which are the most important, because they do not 
involve the sine of the obliquity. 

I shall therefor^ evaluate the tangential stresses when the obliquity is zero. 

The complete potential due both to the moon and to the diurnal rotation is 

r^S = \r'^{n^ + t) (f — cos^ B) + ^^h sin^ ^ cos 2 (^ — • cw^ — e) , 
and the complete expression for the surface of the spheroid is given by 



<j 



g"'=:i(#4-T)(-|— cos^ #)+ii" COS 2€ sin^ 6 cos 2((^— co^). 



jzLence 



a 



S — g-=lTsin 2€sin^ 6 sin 2(^— co^). 

Cv 



Then neglecting r^ compared with m^, and omitting the terms which were previously 
considered as giving rise to secular distortion, we find 






2^ H/ S— g- )=tOTV|™ sin 2€ sin $ cos ^(-|— cos"^ 6) sin 2(<j&— co^), 



wa^- --r~-^^( S—g- )==te?aV|'-^ sin 2€ sin ^(-1— cos^ 6) cos 2{(f>--o)t). 

The former gives the tangential stress along, and the latter perpendiculal" to, the 

meiidian. 

If we put e=|^, the ellipticity of the spheroid, we see that the intensity of the 

S 
tangential stresses is estimated by the quantity wa^,re sin 2€. But we must now find 

a standard of comparison, in order to see what height of tide such stresses vfould be 
competent to produce. 

It appears from a comparison of equations (7) and (8) of Section 2 of the paper on 
'' Tides/' that a, surface traction S^ (a surface harmonic) everywhere normal to the sphere 
produces the same state of flow as that caused by a bodily force, whose potential per 

unit volume is \p]%i\ and conversely a potential W,: is mechanically equivalent to a 



surface traction ( - ] W,-. 



'a 

4 B 2 



Now the tides of the first order are those due to an effective potential wrH^'-^\ 
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and hence the surface normal traction which is competent to produce the tides of the 

first order is ivaHS — Q-l which is equal to tTO^^-r sin 2€ sin^ ^ sm 2(^— cu^). Hence 

the intensity of this normal traction is estimated by the quantity wa%r sin 26, and this 
affords a standard of comparison with the quantity iva\e sin 2e, which was the 
estimate of the intensity of the secondary tides. The ratio of the two is 2e, and 
since the ellipticity of the mean spheroid is small^ the secondary tides must be small 
compared with the primary ones. It cannot be asserted that the ratio of the heights 
of the two tides will be 2e, because the secondary tides are of a hi^^her order of har- 
monics than the primary, and because the tangential stresses have not been reduced to 
harmonics and the problem completely worked out. I think it probable that the height 
of the secondary tides would be considerably less than is expressed by the quantity 
2e, but all that we are concerned to know is that they will be negligeable, and this 
is established by the preceding calculations. 

It follows, then, that the precessional and nutational forces will cause no secular 
shifting of the surface with reference to the interior, and therefore cannot cause any 
such geographical deplacement of the poles, as has been sometimes supposed. 



lie The distribution of heat generated by internal friction and secular cooling. 

In the paper on ''Precession" (Section 16) the total amount of heat was found, 
which was generated in the interior of the earth, in the course of its retardation by 
tidal friction. The investigation was founded on the principle that the energy, both 
kinetic and potential, of the moon-earth system, which was lost during any period, 
must reappear as heat in the interior of the earth. This method could of course give 
no indication of the manner and distribution of the generation of heat in the interior. 
Now the distribution of heat must have a very important influence on the way it will 
affect the secular cooling of the earth's mass, and I therefore now propose to investigate 
the subject from a different point of view. 

It will be sufficient for the present purpose if we suppose the obliquity to the 
ecliptic to be zero, and the earth to be tidally distorted by the moon alone. 

It has already been explained in the first section hov/ we may neglect the mutual 
gravitation of a spheroid tidally distorted by an external disturbing potential wr^8, 

if we suppose the disturbing potential to be wrH S— g~ j, where 'r=a+cr is the equation 

to the tidal protuberance. 
It is shown in (4) that 

S — ■ g-~= |t sin 2€ sin^ sin 2 (<^ — w^). 
If we refer the motion to rectangular axes rotating so that the axis of re is the major 
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axis of the tidal spheroid, and that of % is the earth's axis of rotation, and if W be the 
effective disturbing potential estimated per unit volume, we have 



W=t(;WS--g-)=t(;Tsin 2e.cc2/ (^^) 



It was also shown in the paper on ^'Tides'' that the solution of Sir W. Thomson's 
problem of the state of internal strain of an elastic sphere, devoid of gravitation, as 
distorted by a bodily force, of which the potential is expressible as a solid harmonic 
function of the second degree, is identical in form with the solution of the parallel 
problem for a viscous spheroid. 

That solution is as follows : — 



1 



a 



1% 



(4 



a 



2 2 1^.2\ 



wn 



dW 
clx 






\f. 



with symmetrical expressions for fi and y. 



^. d/W 
bmce — ' 



1 ^AV ^cc 

■— ^ W, the solution may be written 



dx\r'' J r^ clx '/^ 



a\ 



38u 



Lvvv 



, y8=&C.^ y=&c. 



Then substituting for W from (20) we have 



%0T 



a=;— -sin 2€[(8a^—6r^)?/+4irY 



lOT 



fi= ^ sin 2e [(Sa^— 5r^)x+4:Xy^'] 



y=^^sm 2eixyz 



38u 



« « 



• • • (21) 



lOT 



Putting K= .^ sin 2€, we have 



19u 



And 



da 

dx 

d^ 

dx 

dy 

dx 



-Kxy, ~=^K[8a^-{x^+ny^+5z^)'], 



dy 

iK[8a2- (1 5a;""+2/H Sz^)], 

dy 



dl3 



d, 

d: 

d/3 



5K.yz 



■^ 



2Kyz, 



dy 



dy 

2Kzx, 



-Kxy,^=^5Kxzy. . . (22) 



dz 

dy 

dz 



2Kxy 



d7. dy dx ^ dz ^ ' dy dx ^ ^ 



=^ See Thomson and Tait's ' Nat. Phil.,' § 834, or " Tides," Section 3, 
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Now if P, Q, R, S, T, U be the stresses across three mutually rectangular planes at 
X, yyZ^ estimated in the usual way, then the work done per unit time on a unit of 
volume situated at x, y, % is 



clx dy 



dz 



\ dz dy) ' ' \dx dz 



\dy dx 



But P= — j9+2i^-p S=i;( x"+;r )> ^^^^ Q' -^^ '^' ^ have symmetrical forms. There- 
fore, substituting in the expression for the work (which will be called -37), and 
remembering that 



dt 



we have 



da d^ .dj^_^ 

ax dy dz 



i,<j 



.1. CbjLLi 

V dt 



r/^.y /W /M^l Ai^ 

i\dx) \dy I ^\dzjS\dz 



dyV /dy dcA^ (da. d^V 

dyj \dx dzj \dy dx j 



Now from (22) 



2 



dx] \dy ) \dz 



= r2ccy=fr*sin^^[l— COS 4(^—0)^)] . . (24) 



and from (23) 



1 



d^ , djY' fdy doc\^ I da d^\^ 
d^z dyJ \dx dz) \dy dx 



^z^x^-^y^) +[8 (a^—x^—y^) — 5z^J 



9r* sin^ cos^ 6'j-{8a^—5r^'—Sr^ sin^ 6f . 



» 3 



« « 



(25) 



Adding (24) and (25) and rearranging the terms 



1 dM 



K^v dt 



= — |r* sin^' cos 4(^-»a)^) + (8#— 5r^)^-~|r^ sin^ 0[32a^-'Q^%26+ sin^ 0)[. 



The first of these terms is periodic, going through its cycle of changes in six lunar 
hours, and therefore the average rate of work, or the average rate of heat generation, 
is given by 



dK 1 A/;t \^ 

— =4- sin 2€) [{8a^-5rY-ir^ sm^^ 0{S2a^--r%26+ sin^ 0)}] . 



. (26) 



It will now be well to show that this formula leads to the same results as those 
given in the paper on ^^ Precession/^ 

In order to find the whole heat generated per unit time throughout the sphere, we 

r r r 7Tr 

must find the integral i-^^r^ gin 0drd0d<p, from r=a to 0, 0=7t to 0, ^=27r to 0, 
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Thomsoit and Tact, 'Kat. Pliil.,' S 670, 
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In a later investigation we shall require a transformation of the expression for — , 

and as it will here facilitate the integration, it will be more convenient to efiect the 
transformation now. 

If Q3, Q4^ be the zonal harmonics of the second and fourth order. 



Now 



(8a^— 5r^)'-^— fr" siii^ ^[32a^— (26+ sin'^ 9)r^^ 
== (8a'-— Sr^)^- r^ASa^—^i^^ -|(32#_ 28'/-) cos^ 6'— f r^ cos* 0' 
=-l-{320a'^— 560aV^+259r*}-|(112a2— 95r3)r2Q3+ifr*Q4 . 



(27) 



The last transformation being found by substituting for cos^ 6 and cos^ 6 in terms of 
Q.2 and Q,^, and rearranging the terms. 

The integrals of Q^ and Q^^ vanish when taken all round the sphere, a.nd 



4:7rc(? 



Ga^ 



i{320aJ^-560a^i^+259r^7^ mi edrd0d(t>=-^^{^Q-^^+^^] = 



u 



where G is the earth^s moment of inertia, and therefore equal to -^^irwa^. 
Hence we have 



3 



f o 



r^ sin 0drd6d(f)=z-i ^^ sin 2e ) Ca^.f X 19 = 77~--(t sin 2e)^C. 



cU 



u\19 



o 



Sv 



-r, , , ^ 19vco \^ 19vo) ,T , biocc'^ co , ^ 

But tan 2e~ = 2.-—- — -, so tliat -----zz:- cot 2e, 



f^f. 



And the whole work done on the sphere per unit time is | sin 4e.Coi. 

Now, as shown in the first part (equation 5), if J^ be the tidal frictional couple 



-3 



a ^g 

Therefore the work done on the sphere per unit time is J^co. 

It is worth mentioning, in passing, that if the integral be taken from ^a to 0, we find 
that '32 of the whole heat is generated within the central eighth of the volume ; and 
by taking the integral from fa to a, we find that one-tenth of the whole heat is generated 
within 500 miles of the surface. 

It remains to show the identity of this remarkably simple result, for the whole work 
done on the sphere, with that used in the paper on '' Precession. ^^ It was there shown 



* 



Todhuntir's * Functions of Laplaci,' &0.5 p. 13 ; or any otlier work on tlie sabject. 
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(Section 16) that if n be the earth^s rotation, i^ the moon^s distance at any time, p the 
ratio of the earth's mass to the moon's^ then the whole energy both potential and 
kinetic of the moon-earth system is 



■Ciw 



5gT 

2z/'r/* 



Now c being the moon's distance initially, since the lunar orbit is supposed to be 



Also 



Therefore 



V 



r"~\/2 







1 






a^^ 



LW 



-)v%i+v)\ n 







sfl 



5 



according to the notation of the paper on " Precession/^ 

1 
In that paper I also put -=871^12^^ 



f^ 



Therefore ^. i=^. 

2v r p 



2 ^oj^o 



And the whole energy of the system is |C w 



I /74? UH') n f'x 

Therefore the rate of loss of energy is— CM^— 4-Ti'^'^o;^ )• 



dt 



dt 



dfb ^ , dP ^ 12 

But -j-zzz-^-^^ and as shown in the first part (19), jii%-37=— , also -i;^=/3. 



^ , ^..^ ^. .... .... ..... j,^.. ,..,, ^..^^^ ^ 

Therefore the rate of loss of energy is j^(fi^— i2) or jS^i which expression agrees 
with that obtained above. The two methods therefore lead to the same result. 

I will now return to the investigation in hand. 

The average throughout the earth of the rate of loss of energy is J^co ~- l^ra^, which 
quantity will be called H. Then 



H= 



.3 



.2 



4:7ra^ 



0) = -~ • 4Ma^ . 4— sin ie.o)= ¥iva^ . — sin 4 e . w. 



M 



tC 



S 



8 



Now 



H ■""- sm 2e a*=-^r- cot 2€. rr sm^ 2€.t(^a"=Tir--l-t^a'".---- sm 4€.a) = 77;, 
u\19 / 6g 19 19 6 g 2^9 



Hence (26) may be written 



dt "^ 19 



8-61 



r 



a 



9>1 S 



3 
2 






,C(^. 



sin^^^32-(26+sin^^) 



■''^AS^ 



,^. 



(28) 
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This expression gives the rate of generation of heat at any point in terms of the 
average rate, and if we equate it to a constant we get the equation to the family of 
surfaces of equal heat-generation. 

We may observe that the heat generated at the centre is S-^-g times the avera^ge, at 

of the average. 



1 1 

the pole — - of the average, and at the equator -^ 



n 



The accompanying figure exhibits the curves of equal heat-generation ; the dotted 
line shows that of ^ of the average, and the others those of ^, 1, 1|-, 2, 2|, and 3 times 
the average. It is thus obvious from inspection of the figure that by far the largest 
part of the heat is generated in tlie central regions. 



Eio^ 2. 




The next point to consider is the effect which the generation of heat will have on 
uiiderground temperature, and how far it may modify the investigation of the secular 
cooling of the earth. 

It has already been shown''' that the total amount of heat which might be generated 
is very large, and my impression was that it might, to a great extent, explain the 
increase of temperature underground, until a conversation with Sir W. Thomson led 
me to undertake the following calculations : — 

We will first calculate in what length of time the earth is losing by cooling an 
amount of energy equal to its present kinetic energy of rotation. 

The earth's conductivity may be taken as about '004 according to the results given 
in Everett's illustrations of the centimeter-gram -second system of units, and the 
temperature gradient at the surface as 1^ C. in 27-| meters, which is the same as 
1^ Fahr. in 50 feet — the rate used by Sir W. Thomson in his paper on the cooling of 

the earth, t 

4 
This temperature gradient is - — r-^ degrees C. per centimeter, and since there are 

31,557,000 seconds in a year, therefore in centimeter- gram -second units. 
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The heat lost by] 4 4 

aa-^i-h T^^■,. o,.v..,,.. f = ^aith's siiiface in square centimeters x 777 x -r~~z~Zo x 31557 x 10^ 
udriin pel annum j ^ lU ii x lU 

= earth's surface x 45*9 (centimeter-gram-second heat units). 

Now if J be Joule's equivalent 

Earth's kinetic eneroyi /s o i^r / o \ 

of rotation in heatl=-| — ^=— (f)'^( f—— ), where C = fMci^ 

^J J \ ^ 



units 



earth's surface x -^r~-{f)\, where eQ=f-^ =^1^ 

3J g 

^, , . (5'6)x (6-37)'- X 1016 X (-4)3 for a=:6-37 x 10^ centimeters, 
earth s surface x — ^ £„„_^,-_„™_. — ^ — l., 

3 X 4-34 X 10* X 232 J = 4-34 x lO'^gramcentim. 



= earth's surface x 1'2 x 10^^ nearly. 



and 10 = 6^. 



Therefore at the present rate of loss the earth is losing energy by cooling equivalent 

1*2 X 10^^ 
to its kinetic energy of rotation in -~^^^—~ — =262 million years. 

If we had taken the earth as heterogeneous and C=^Ma^ we should have found 218 
million years. 

We will next find how much energy is lost to t lie moon- earth system in the series 
of changes investigated in the paper on '^ Precession." 

In that paper (Section 1 6) it was shown that the whole energy of the system is 

^MaHn^—-T- -;], where v is earth -r- moon, r moon's distancCj n earth's diurnal rotation. 



Hence the loss of energy =i-Ma%o^ 
to tIq, and r from r to Tq. 



[_\nj 2vn^ \r r^/j 



, while n passes from n 



Now — ^=--- — ^ a=-~~---~:^"a=8'84a, taking z^=82, and p-^=232. 
2pnQ^ Sv \5?'?.o%/ 32 X 82 ^ 5^o^ 

n I) 
If D be the length of the day, — =-^ ; and if U be the moons distance in earth's 



radii, then 



loss of energy = 



/ir\ JlJ 



T) \2 / 1 1 

^10) — 1_8-84(- — 



X earth's present he. of rotation, 



But in the paper on ^^ Precession'' we showed the system passing from a day of 5 hours 
40 minutes/"" and a lunar distance of 2*547 earth's radii, to a day of 24 hours, and a 
lunar distance of 60*4 earth's radii. 

* A recalculatioii in the paper on " Precession " gave 5 hours 36 minutes, but I Lave not thought it worth 
while to alter this calculatioD. 
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Now 24-T-5|=4-23, and (2-547)""i — (60-4)-"^=-376. 

Therefore the loss of energy=[(4'23)^—l— '376 X 8*84] X earth's present k,e. 

= 13*57 X earth's present he. of rotation. 

Hence the whole heat, generated in the earth from first to last, gives a supply of 
heat, at the present rate of loss, for 13*6x262 miUion years, or 3,560 million years. 

This amount of heat is certainly prodigious, and I found it hard to believe that it 
should not largely affect the underground temperature. But Sir W. Thomson pointed 
out to me that the distribution of its generation would probably be such as not 
materially to aJfect the temperature gradient at the earth's surface ; this remarkable 
prevision on his part has been confirmed by the results of the following problem, which 
I thought might be taken to roughly represent the state of the case. 

Conceive an infinite slab of rock of thickness 2a (or 8,000 miles) being part of 
an infinite mass of rock ; suppose that in a unit of volume, distant x from the 
medial plane, there is generated, per unit time, a quantity of heat equal to 
5[320a^^— 560aV+259cc^]; suppose that initially the slab and the whole mass of rock 
have a uniform temperature V; let the heat begin to be generated according to the 
above law, and suppose that the two faces of the slab are for ever maintained at the 
constant temperature V ; then it is required to find the distribution of temperature 
within the slab after any time. 

This problem roughly represents the true problem to be considered, because if we 
replace x by the radius vector r, we have the average distribution of internal heat- 
generation due to friction ; also the maintenance of the faces of the slab at a constant 
temperature represents the rapid cooling of the earth's surface, as explained by Sii 
W. Thomson in his investigation. 

Let S- be temperature, y thermal capacity, k conductivity ; then the equation of heat- 
fl.ow is 

/|z.^f|+|)r320a^-560aV+259irn 
Let 320^=2L, 560^=:12M, 259 ^^30N, and let 'the thermometric conductivity 



/c=-~. Then 

y 

clS^__ d? 

dt dx- 



'^+La^c^^— MaV+Ncc®— R" 



Let the constant Rr3:(L— M+N)a®, and put 






Then when ir=±a, t//=5^. 
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Since Lj M, N, R are constants as regards the time, 

— -^=.K — ~, 

(It ■ dx^ ' 

i//=y ■— SPe"""^^'^ COS qx is obviously a solution of this equation. 

Now we wish to make .9^::=V5 when. ^=:--^a^ for all values of t; since i//=:5- when 

xi^ia, this condition is clearly satisfied by making g=:(2?'+l)5". 

Hence the solution may be Avritten, 

^=:V-[Lc^i7;2-M:aV+Nx«-R]-SP2;+ie-^'t-1f'^ cos(2i+l)-^ . (29) 

and it satisfies all the conditions except that, initially, v/hen ?^=:0, the temperature 
everywhere should be V. This last condition is satisfied if 

SP2/+1 COS (24+l)J-=R-La%HMaV-N«« 

for all values between (r=:i:a. 

The expression on the right must therefore be expanded by Fourier's Theorem ; 
but we need only consider the range from x^:=:a to 0, because the rest, from x=-0 to 
— a, will follow of its own accord. 



m^- 



Let Y^ ^^''' ; let ot be written for - : let M'=— , N' = -, and R^^R^, 
Then 



1j5' 



and this has to be equal to :SP2h-iCOS (2{+1)x fron^ X^o" ^^ ^' 



'^ 



Since 



TT 



and 



2 cos (2^*+l)xcos (2y+l-)x^^X"^ un.lessy=:?*3 



2 COS ^(2{+ l)x^7x^^ i'^^i^* 







Therefore 



«/^ ^"- 



ir^P,,, , = ^ '^ [K ~ Lx" + M'x* - N'x'l cos (2i + 1 )x'^X- 



■US" J 







Now 
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1 



2i-\-l 






dv^ 



X^^-sin (2^+l)x+-^ cos (27:+l)x--^-jsin (2i+l)x 



21 + 1 






(2i + l) 



(2^ + 1) 



- COS (2i+ l)x+&c. 



7r 
2 







(-y 



1- 






d 



i 



+._^^_&c. 



2* + l_' (2^ + 1)2 (2* + !)'' 



2? 



If therefore /*(;)() be a function of ;)( involvixig only even powers of ;>(, 



f/(x)cos(2^+l)x^X^ 



(-y 



(2^ + 1) 



1 + 



1 dV 



2^^~l <:fe 



/(^). 



This theorem will make the calculation of the coefficients very easy, for we have 
at once 



m 



2a^ ■^'-'■+'^ 



2i + l 



|R'„Lt^2_}.mV'-NW' 



1 



/a- , ,,-r— 2L+4.3MV-6.5N'ot*1 
+ /o-ViMr4.3.2.1.M'-6.5.4.8N'w21 

- TTTT^J --■6.5.4.3.2.lNn]. 
(2^ + 1)6'- Jj 



Substituting for R', L, M', N^ their vahies in terms of - we find 



-L 2?+: 



(2^ + l)V Z; 



19 



1988 



"0 + 



6216 



(2i + iym'^ ' (2^^-l)Vl 



Then putting for zu its value, viz.: ^ of 3*14159, and putting i successively equal to 
0, 1, 2, it will be found that 



W 



Pi==f-(120-907), P3-f (1-107), P5 



17 



(•048y 



So that the Fourieu expansion is 



120'907 cos — +1-107 cos-:~-048 cos 



9, 



a 



2a 



2a' 



which will be found to differ by not so much as one per cent, from the function 
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m ^-{tf)-M' -(Sf)+'»{'-(Sf '• 



to which it should be equal. 

Then by substitution in (29) we have as the complete solution of the problem 
satisfying all the conditions 

3=Y+~i 1-6 ^"'^^ 120-907 cos- 
h i\ J 2a 

+ (l--e ^-'^ J 1-107 cos -^ il—e ^^^ J '048 cos — 

The only quantity, which it is of interest to determine, is the temperature gradient 

at the surface, which is equal to —-7- when x=:4zCL 
Now when x= J^a, 

" I ^ is a small fraction, we have approximately 
^^^|a:^^^_7^YJl20-907-9X3-321-25X-240l; 



dx h 2 \2a, 



tc 1 

and since 7=- 



fAfjW')- 



K\ 



This formula will give the temperature gradient at the surface when a proper value 
is assigned to ^, and if t be not taken too large. 

With respect to the value of t, Sir W. Thomson took k=400 in British units, the 
year being the unit of time ; and a=21x 10^ feet. 

X.JLtJxl'w't? 

2^j =4^ l^\2T^a^TJ =ro^. nearly, 

"I =rT^; if therefore t be 10^ years, this fraction is 2^-. Therefore the 

solution given above will hold provided the time t does not exceed 1,000 million years. 
We next have to consider what is the proper value to assign to J. 
By (27) and (28) it appears that Ja* is y-i-^j- of the average heat generated 
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throughout the whole earth, which we called H. Suppose that p times the present 
kinetic energy of the earth's rotation is destroyed by friction in a time T, and suppose 
the generation of heat to be uniform in time, then the average heat generated through- 
out the whole earth per unit time is 



Therefore 



~~=: . ¥Kahir? -^ earth^s volume. 



5JT ' T""^^^ 2~5 j^ waeQ, 



m in 

Where e^ is the ellipticity of figure of the homogeneous earth and is equal to f -—, 

which I take as equal to 2-1^2^. 
Hence 

and 

dj% 16 X 85 /ttX^ w p6q t 

But y=zsw^ where s is specific heat. 
Therefore 

d^ 1707T^ pe^ 1 t 
'^dx~ 9500"~s"jT 

The dimensions of J are those of work (in gravitation units) per mass and per scale 
of temperature, that is to say, length per scale of temperature; p^ 6q, and 5 have no 
dimensions, and therefore this expression is of proper dimensions. 

Now suppose the solution to run for the whole time embraced by the changes 
considered in ^^ Precession," then t=T, and as we have shown |7=rl 3*57. Suppose 
the specific heat to be that of iron, viz.: |-. Then if we take J = 772, so that the 
result will be given in degrees Fahrenheit per foot, we have 



clx 


■ 950 


13-57 
^232x 


x9 

772 


1 









2650* 

That is to say, at the end of the changes the temperature gradient would be 1^ Fahr. 
per 2,650 feet, provided the whole operation did not take more than 1,000 million 
years. 

It might, however^ be thought that if the tidal friction were to operate very slowly. 
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so that the whole series of changes from the day of 5 hours 36 minutes to that of 
24 hours occupied much more than 1,000 miUion years, then the large amount of heat 
which is generated deep down would have time to leak out, so that finally the 
temperature gradient would be steeper than that just found. But this is not 
the case. 

Consider only the first, and by far the most important, term of the expression for 
the temperature gradient. It has the form J (1 — e""^^), when i5::^T at the end of the 

series of changes. Now ^ varies as T ^, and ^^ has its maximum value imity when 

T=0. Hence, however slowly the tidal friction operates, the temperature gradient can 
never be greater than if the heat were all generated instantaneously ; but the tem- 
perature gradient at the end of the changes is not sensibly less than it would be if all 
the heat were generated instantaneously, provided the series of changes do not occupy 
more than 1,000 million years. 

III. The forced oscillations of viscous, fluid, and elaMic spheroids. 

In investigating the tides of a viscous spheroid, the effects of inertia were neglected, 
and it was shown that the neglect could not have an important influence on the 
results.'"' I shall here obtain an approximate solution of the problem including the 
effects of inertia; that solution will easily lead to a parallel one for the case of an 
elastic sphere, and a comparison with the forced oscillations of a fluid spheroid will 
prove instructive as to the nature of the approximation. 

If W be the potential of the impressed forces, estimated per unit volume of the 
viscous body, then (with the same notation as before) the equations of flow are 



dx dx \dt dx dy ^ dz) 

-?+&c. = 0-?+&c. =0 > 

dy dz 

da. d/3 j_ ^7_ /^ 

dx dy dz _j 



(30) 



The terms — i^f -_,+ &c. j are those due to inertia, which were neglected in the paper 

on ^^ Tides. '^ 

It will be supposed that the tidal motion is steady, and that W consists of a series of 
solid harmonics each multiplied by a simple time harmonic, also that W includes not 
only the potential of the external tide-generating body, but also the effective potential 
due to gravitation, as explained in the first part of this paper. 



# u m 



Tides," Section 10, 
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The tidal disturbance is supposed to be sufficiently slow to enable us to obtain a first 
approximation by the neglect of the inertia terms. 

In proceeding to the second approximation, the inertia terms depending on the squares 

and products of the velocities, that is to say, '^(^T;+i87 +7T")^ ^^J ^^ neglected com- 

\ U,ih CI If Cviy / 



clot 



pared with t(;~. A typical case will be considered in which W=Ycos (vt-\~e), where 

CvZ 

Y is a solid harmonic of the i^^' degree, and the e will be omitted throughout the 
analysis for brevity. Then if we write 1= 2(^+1 )^+l> the first approximation, when 
the inertia terms are neglected, is 



a= 



Iv 



<^+^ )a^__ (^+;^j(^^+^) ;^^' 



L2(^~1) 



2(2i+l) 



dx 



% 



i.2i+3 



2^+l' 



dx 



{r-2'-iY) 



cos vt'" 



(31) 



Hence for the second approximation we must put 



doL wv 
clt Iv 



sm vt. 



And the equations to be solved are 



dp o . dY^ , . iDv 

dx If \\Ji{% — \) 



f(^ + 2) o (^■fl)(2i + 3) , 



a 



2 



2(2^ + l) 



dx 



*^ 



^ 



2^ + 1' 



.,2i+3. 



d 



dx 



(r-3^-iY) 



sin vt 



= Y ■ (32) 



dp 
dy 



+ &c. =0, 



d^ 

dz 



4- &c. =0 



These equations are to be satisfied throughout a sphere subject to no surface stress. 
It will be observed that in the term due directly to the impressed forces, we write Y' 
instead of Y ; this is because the effective potential due to gravitation will be different 
in the second approximation from what it was in the first, on account of the different 
form which must now be attributed to the tidal protuberance. 

The problem is now reduced to one strictly analogous to that solved in the paper on 

doc 
" Tides;" for we may suppose that the terms introduced by ivy, &c., are components of 

bodily force acting on the viscous spheroid, and that inertia is neglected. 

The equations being linear, we consider the effects of the several terms separately, 
and indicate the partial values of a, /3, y, p by suffixes and accents. 

First, then, we have 

clY' 



— •~-f-^+7;V^an+-T~ cost'^^0, &c., &c. 

dx ^ dx 



*• '* Tides," Section 3, equation (8), or Thomson and Tait, ' Nat. Phil.,' § 834 (8). 
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The solution of this has the same form as in the first approximation, viz. : equation 
(31), with a^ written for a, and Y' for Y. 

We shall have occasion hereafter to use the velocity of flow resolved along the radius 
vector, which may be called p. Then 



X 



P0=«'0 



y 



+^o~ +ro:; 



Hence 



^0^ 



1_ 
Iv 



"2(^-1) 



)-— cos vt 



« » 



(33) 



Then observing that Y'-r-r' is independent of r, we have as the surface value 



Secondly, 



a^+i i(2i-hl) Y' 

Po — T~ T^p—rz'-r QO^ vt , 
' ^ iv 2(^ — 1) r* 



_% +,V Vo+^'-i±|^ ? sin vt = 0, &c., &c. 



• •«»«t0d 



dx 



Iv 2(4-1) dx 



(34) 



(35) 



This, again, may clearly be solved in the same way, and we have 



^0 = 



lova^ i(i-{-2) 



X^+2) , (^+i)(2^+3) ; 



[2(i-l) 



2 



(2^ + 1) 



clY 



I 



dx 2i + l 



7 



.2i+3l(Y7."-^J~l) I gill ^^ (36) 

CvX J 



and 



f^^ivva^ i{i-\-2) 
^0"~lv •2(^-1) 



r /•■) 



#(^4-2)^6--^(^^--l)r- 
2(^-l)"■ 



Y 

- sin 'y^ 

7' 



and its surface value is 



p^ z=L wva 



_ ^^(^ + 2)(2i+l) Y^, 



Thirdly, let 



[2Iu(i-l)P r 



: sm vt 



u 



loj; 



Y 



Iv 2{2i+l) 



sm ■vi 



(37) 



(38) 



(39) 



So that U is a solid harmonic of the i^^' degree multiplied by a simple time harmonic. 
Then the rest of the terms to be satisfied are given in the following equations :— 



dp 

dx 

dp 



-{-VV OL — 



(i+l)(2.-+3),-^f+2»-«">r--)' 



di/ 



+ &c. =0, 



dp 

dz 



dx 



+ &c. =0 



dx 



= 



1 



V 



/' • • « 



^ 



(40) 



These eqiiations have to be satisfied throughout a sphere subject to no surface 
stresses. The procedure will be exactly that explained in Part I., viz.: put a=:a'+a^, 
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jS=/3'+Aj T='y'+7/?P=p'+P/? ^^^^^ fi^^ ^^ ^'j t'?P' ^W functions which satisfy the 
equations (40) throughout the sphere, 

Diflferentiate the three equations (40) by x, y, z respectively and add them together, 
and notice that 



(i+l)(2^+8) 



and that 



d I AV>\ , d 
dd\ dx I dy 






d 
dz^ 



+2i^|/r^^+4Ur-2^-i)) + ~ 



(ajJU 



dx 



dai! d^' dj^ 

•M^MM* MiMaNaHB (^•■■^■n \ / 

dx dy dz ' 



then we have v^2^'=^0^ of which j9'=0 is a sokition. 
Now if V^^ be a solid harmonic of degree % 



V h^''Yu=m{27i+m+ 1)^^'"^,, 



Hence 



r 



a' 



dU_ 

dx 



2. 



dV 



1 



4(2i + 3) dx 



Si+^ 



d 



dx 



i(XJr-.^-i)=^2 



,.2i+5 g 



!> . 



2(2^ + 5) dx 



{JJr-^i-i) 



ma"^ 






d 
dz' 



0, 



t • • 



(41) 



Substituting from (41) in the equations of motion (40), and putting |)'=0, our equa- 
tions become 



3 



V^iVOL — 



, (^+1) M 



2'^ + 5 dy^ '\ 



4 dx 



^» t « • • ItXjWl 



of which a solution is obviously 



1 fi^ + l .<iU 



V 



-r 



4. 



% 



d 



^4 6^0^ 2i + 5 

/3'==&C., y' = &C. 



4_ _Jl__«,2e+5 J:l/TT^-2i-l\ 



<ia? 



> . 



. . (43) 



It may easily be shown that these values satisfy the equation of continuity, and 
thus together with p' = they are the required values of a, fi\ y\ p\ which satisfy 
the equations throughout the sphere. 

The next step is to find the surface stresses to which these values give rise. The 
formulas (13) of Part I. are applicable 



vt^ =v{ax+ ^'y +yz) 

— ^'(^ + 1) 4 jj _ H'^ + 1) 4TT , 

'^ 4 2t; + 5 

4 D 2 



4(2^ + 5) 
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rii 



Then remembering^ that 

o 



XU: 



1 



dV 



sz:. 



2i+l dx 



•2^+3£(5,-2i-lU) I 



ml 
nave 



V 



d^\. J{i-\-l){2i^l) 

dx~^ ' 4(2^+5) 



,.4 



dJ] . 4r3 .fZU 



cm 



.3: 



45«3 



d 



dx 2i+l' 






cZa; 



cw 



f& 



• • « • ♦ • 



(44) 



Again, by the properties of homogeneous functions, 



d 



v\ r-~-\y=v[xj^+yj+Zjy-va 



dr ""/"' ^\ dx "^ dy dzj 

=^ ^i --^-r^''!" +-^T^ . . . . (45) 

4 fZ^ 2^ + 5 dxr ^ ^ 

Alsoj9^=0. 

Then adding (44) and (45) together, we have for the component of stress parallel to 
the axis of ^ across any of the concentric spherical surfaces^ 



'Fr='-'p'x-\-v 



(i + l)\JU 



r-T-— 1 06 4 



dr 



dx 



by (L3)j Part I. 



d 



2 



'^+2iT^^"'^^^^^^ ^y (^^) ^^^ (^^)- 



And at the surface of the sphere^ where r=a, 



J^--lllll£^l+2 



-1+1 



dJI 

dx 



% 



"^2^ + 5^ 



i+2 



d 



rf. + %lLL.-9,i'-l'\J\ 



Civy 



* e * « 



(46) 



The quantities in square brackets are independent of r^ and are surface harmonics of 
orders i—l and i-\-l respectively. 

Let 



Where 



-^i- 



«-i 



(^ + 1) 



3 



a 



e+2 



F: 



...cZU" 



A. —A. 



r 



-i+i 



dx 



, A,'+i= 



^ 



2i + 5 



a 



2 -J- 2 



d 



^J-^2JZ./rp-2i-lJJ\ 



dx 



-" J 



■ (47) 



Also let the other two components G and H of the surface stress due to a', jS^ y , p' 
be given by 

Then by symmetry it is clear that the B^s and C's only differ from the A's in having 
y and % in place of x, 
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We now have got in (43) values of ot', j^\ y\ which satisfy the equations (40) 
throughout the sphere, together with the surface stresses in (47) to which they 
correspond. Thus (43) would be the solution of the problem, if the surface of the 
sphere were subject to the surface stresses (47). It only remains to find a^, ^^, y^, to 
satisfy the equations 



_^|^+.vV=o, -%+&c.=0, -^^' 



dx 



dy 



dz 



f&c. = 



(48) 



throughout the sphere, which is not under the influence of bodily force, but is subject 
to surface stresses of which A^-^i+A.+i, B^-^^+Bi+i, C!i«i+Q-i-i ^^^ the components. 

The sum of the solution of these equations and of the solutions (43) will clearly be 
the complete solution ; for (43) satisfies the condition as to the bodily force in (40), 
and the two sets of surface actions will annul one another, leaving no surface action. 

For the required solutions of (48), Sir W. Thomson's solution given in (15) and (16) 
of Part I. is at once applicable. 

We have first to find the auxiliary functions ^/_o, *i corresponding to A^^i, B/_i, 
Q„i, and "^i, #^+3 corresponding to A^+i, B^+i, Q+p It is easy to show that 



^. .=0, 4)^+3=0, 



and 



-5/ 



=:-"a'*+^ 



% 



2^ + 5 



dx 



. d .__ 1 

/v»^? + 0____ //l',""j<.2"*"l I j j L I 
l'"A/ J 



£ 

dy 






dz 



xi 



2^ + 5 



$.:= 



-.2?"+ 1^2 + 3 



(i + l)2 



dx\ dx I d^f 



+ 



dz' 






We have next to substitute these values of the auxiliary functions in Thomson's 
solution (15), Part I. It will be simpler to perform the substitutions piece-meal, and 
to indicate the various parts which go to make up the complete value of a^ by accents 
to that symbol. 

First For the terms in a^ depending on A,-^i, %^2^ ^i, we have 



1 



a = 



d^. 

z 






j—i 



va'-^[2(i-2){i-l)(2i-T) dx ' i-2 



a-' 



i(i + iy dJJ (i+iy dJJ 



V [4:{i-l){i-2) dx 2(i-2) dx 

«^- (^+1)^ dJJ 

V 4(i—V) dx ' ' ' ' 



' (^9) 



(Note that i—2 divides out, so that the solution is still applicable when i=2). 
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Second. In finding the terms dependent on A^^^, '^i^ ^i^^ it will be better to 
subdivide the process further. 

^ ' va'2 1^ ^ dm 

^v 2I(2^ + 5) ^^ ^^& iou; 



(ii) a/^=: ^^ { -^±1^ + 1a,,i7-- 



+1 



v[ I(2t4-5) ffe^ ^ 2^ + 5 ffo"^ ^J 



Then since 



therefore 



(^+3)(^+l)■~I=^H4^+3-2^^-4^-3 = -^^ 



a'"=_'i^^_l_-,.2.-+s5i(,,-a--]XJ) . ...... (51) 



/ 



V I(2i + 6) rfe 



This completes the solution for a^. 

Collecting results from (49), (50), and (51), we have 



a =06, +oe +06 






V 



[4(^-^lf rfe 21(2^+5) ^ ^cfe^I(2^ + 5) dx^ ^^^ ' ^ ^^ 



Then collecting results, the complete value of a as the solution of the second 
approximation is 

06:=:06Q-|~otQ -j-a -j-Ot . 

So that it is only necessary to collect the results of equations (31), (with Y' written 
for Y), (36), (43), and (52), and to substitute for U its value from (39) in order to 
obtain the solution required. The values of ^ and y may then at once be written 
down by symmetry. The expressions are naturally very long, and I shall not write 
them down in the general case. 

The radial velocity p is however an important expression, because it alone is 
necessary to enable us to obtain the second approximation to the form of the spheroid, 
and accordingly I will give it. 

It may be collected from (33), (37), and by forming p and p^ from (43) and (52). 

I find then after some rather tedious analysis, which I did in order to verify my 
solution, that as far as concerns the inertia terms alone 
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! - sin vtmr^-Ma^r^-^ata'^] , 



Where 



®: 



i(i + 1) 
2.4(2i + 5)I' 

2 



21 



i3_!!(!±lK2!!±10i±?). ^^ 
*'— 4(^-1) (2i + 5)F ' 



.(i + 2\^ ('i + l)(2i + 3) 



(2^+1) (2^ + 5) J 2.4(^-1) (2^-f 1)1 



If Qt be reduced to the form of a single fraction, I think it probable that the 
numerator would be divisible by 2^+lJ but I do not think that the quotient would 
divide into factors, and therefore I leave it as it stands. 

In the case where ^= 2 this formula becomes 

p=~ ~ sin ^'«^-=f^{19r*-I48aV+287a'^■}, 



y2 Q. 



213.192 



which agrees (as will appear presently) with the same result obtained in a diflferent 
way. 

I shall now go on to the special case where ^=2, which will be required in the 
tidal problem. 

From (39) we have 



U 



tov 



1 



V 2.5.19 



rY sm vt. 



From (36) 



iDva'- 



S = 



V 



2 



192 



^ , 3.7 MY 2 ^d 
2.5 I dx 5 dx 



r^ ] "z rr'-7- (Yr ^) 



sin vt. 



From (43) 



From (52) 



wv 



1 



a 



v^ 213.5.19 



- a(-V JL , ^.4fc f\ (.V /T7- . ^\ 



dx 



dx 



sin vt. 



'lO'VCt' 



a 



v^ 2'13.5.19' 



(5. 9 7aH V'. 4r^)£- Y'^^-r-^Y) 



Adding these expressions together, and adding a^, we get 



, wv 

a=aoH — ¥'^ 



1 



v'- 2^.3.5.192 



dY 



(5.287a4'-37.4.raV^+9.19'r^^)-- 



C-& 



8 



d 



(2.37a^-"19r^)?^^— (Y^^""^) 



(^03 



-J 



sm vt 



(53) 



and symmetrical expressions for ^ and y. 



-2! 



In order to obtain the radial flow we multiply a hj -, /3 by -, y by -^, and add, 



and find 



ni^ 1 Y^ 

/'=Po+-;;r-^iJi^(287a*-4.37aV2+i9r*)-sin(^«+e) . . . (54) 



the e which was omitted in the trigonometrical term being now replaced. 
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The surface value of p when r=.a is 

p=p^^—^^^^^sm{vt+e) ....... (55) 

where p^ is given by (84). 

If we write -—^Tr—e for e we see that a term Y sin (vt—e) in the effective disturbing 

potential will give us 

ivvcc^ 79 Y 
p=p^-^-^^^~j~~^^G08{vt-e) . (56) 

Now suppose tt?r^S cos t^<^ to be an external disturbing potential per unit volume of 
the earth, not including the effective potential due to gravitation, and let r=a+o-^ 
be the first approximation to the form of the tidal spheroid. Then by the theory of 
tides as previously developed (see equation (15), Section 5, '^ Tides") 

-'= - COS € cos ivt ■— e), where tan e= 7; — ■ 

C6 g . ^ ^ 2gmD 

Then when the sphere is deemed free of gravitation the effective disturbing 
potential is wrH S cos vt'—^~\; this is equal to --iDr'^mi e S sin (vt'—e). 

Then in proceeding to a second approximation we must put in equation (56) 

Y=— tyrosine S. 

Thus we get from (56)^ at -the surface where r=a^ 

p^pQ-i ^ • ,— — 2 sm € b COS (1;^ — e) ...... (57) 

To find Pq we must put r=a-^(T as the equation to the second approximation. 
Then Pq is the surface radial velocity due directly to the external disturbing potential 
tvr^B cos vt and to the effective gravitation potential. The sum of these two gives an 

effective potential tvrHS cos vt—Q^\ which is the Y' cos vt of (34). 



Then pQ is found by writing this expression in place of Y^ cos vt in equation (34), and 
we have 






Substituting in (57) we have 



-^( S COS ««-S- + -j^- ^j^, sm e S cos [vt-e) ) . . . . (58) 



With the tides of a viscous spheroid. 575 

Then since tan e=^— — , therefore '-—-"='" cot e. and (58) becomes 

p=a- cot e( S cos vt — ^- +77^^— cos e S cos (vt-^e) . 
'^ g \ ^a ' 150 g ^ V 

But the radial surface velocity is equal to jj, and therefore j7^=p} so that 



"17 -{-V cot e.or=^a - cot e( S cos t^i^ + .7- - cos e S cos {vt — e) ) . . . (59) 

Then if we divide cr into two parts, cr\ a\ to satisfy the two terms on the right 
respectively, we have 

cr S / v^ 

- = COS €•- COS (vt — e), 
a. g ^ ^' 

which is the first approximation over again^ and 

G-^^ S 79 ^2 

—==: cos e •" -7-37, — cos e cos (vt — 'Ze), 

a g 150 g ^ ^ 



Therefore 



a 



s r ^9 v^ 1 

coS€---j cos (^'i^— e) + — -" COS ecos (^'^~-2e) V .... (60) 



This gives the second approximation to the form of the tidal spheroid. We see that 
the inertia generates a second small tide which lags twice as much as the primary one. 

Although this expression is more nearly correct than subsequent ones, it will be well 
to group both these tides together and to obtain a single expression for cr. 

Let 



tan;X= 




tA— sin € cos e 






Then 



l+-i¥o^cos"-e 



a g cos ;j^\ ' -^ ^ " g ^ ^ ^' ^ ' 



COS G I V?^ 

This shows that the tide lags by (e+x)' ^'^^ ^^ ^^ height ( 1 +1-5-0-; cos^^ e ) of the 



equilibrium tide of a perfectly fluid spheroid. 






By the method employed it is postulated that 1-5-0 ~ is a small fraction, because the 
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effects of inertia are supposed to be small. Hence x niust be a small angle, and there 
will not be much error in putting 

^=1-^^^— sin e cos e, and sec ;)(= 1. 
S 

Then we have for the lag of the tide (e+i-s^^sin e cos e), and for its height 
cos eM +^5%^ cos^ e j. 
Let 7) be the lag, then 

79 ^ 



^ = e+f5o-r Sine cose, 

S 



whence 



^3 



Also 



e=7] — i-g-Q— sin 7) cos t) very nearly, 

S 



cos e=^ cos 7][ 1+1-5-0"- sin^ t] 



and 



/y3 \ / q-)^'' 

COS e( 1 + 1-5-0- cos^ e = cos 7j 1 +1-5-0 ;; )• 



Hence (61) becomes 



O- S / ^ . „r. v^ 



wher 



e 



,^ = g cos 7j(^l+-f5%~ ) cos {vt-rj), 



Tj— -i-g-Q— sin 7j cos 7j = arc tan 



>..,.,, (62) 



^gawj 



-/ 



This is probably the simplest form in which the result of the second approximation 
may be stated. 

From it we see that with a given lag, the height of tide is a little greater than in 
the theories used in the two previous papers ; and that for a given frequency of tide 
the lag is a little greater than was supposed. 



The whole investigation of the precession of the viscous spheroid was based on the 
approximate theory of tides, when inertia is neglected. It will be well, therefore, to 
examine how far the present results will modify the conclusions there arrived at. It 
would, however, occupy too much space to recapitulate the methods employed, and 
therefore the following discussion will only be intelligible, when read in conjunction 
with that paper. 

The couples on the earth, caused by the attraction of the disturbing bodies on the 
tidal protuberance, were found to be expressible by the sum of a number of terms, 
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each of which corresponded to one of the constituent simple harmonic tides. Each 
such term involved two factors, one of which was the height of the tide, and the other 
the sine of the lag. Now if e be the lag and v the speed of the tide, it was found in 
the first approximation that tan €=19u^'-^2^a^(;, and that the height of tide was pro- 
portional to cos € ; hence each term had a factor sin 2e. 

But from the present investigation it appears that, with the same value of e, the 

height of tide is really proportional to cos €(l + i~5-o~" cos^ e j ; whilst the lag is 

€+1-5-0"" sin € cos €, so that its sine is ( 1 + 1-5-0 - cos'^ e ) sin e. 

Hence in place of sin 2e, we ought to have put sin 2e( l + i~5^~cos^ ej , or 

sin2efl+^f-cos2e 



di dN dP 
Thus every term in the expressions for ^. -77; 77 should be augmented, each in a 

Ci/v CvJ/ U/u 

proportion depending on the speed and lag of the tide from which it takes its origin. 

''' In the paper on '' Precession," two numerical integrations were given of the 
differential equations for the secular changes in the variables ; in the first of these, 
in Section 15, the viscosity was not supposed to be small, and was constant, in the 
second, in Section 17, it was merely supposed that the alteration of phase of each tide 
was small, and the viscosity was left indeterminate. It is not proposed to determine 
directly the correction to the first solution. 

The correcting factor for the expression sin 2e is greatest when e is small, because 
cos^ e may then be replaced in it by unity ; hence the correction in the second 
integration will necessarily be larger than in the first, and a superior limit to the 
correction to the first integration may be found. 

We have tides of the seven speeds 2{n—n), 2n, 2{n-\~n), n—lfL, n, n+2/2, 2/2; 
hence if the viscosity be small, the correcting factors for the expressions sin 46^, sin 4e, 

sin 4 £2, sin 2e'i, sin 2e', sin 2e'2, sin 4e'' are respectively l+fl;, multiplied by the 
squares of the above seven speeds. 

Then if X^—, the seven factors may be written 



n^ 



•^ 



1+W~" multiplied by {l—Xf, I, (1+^)^ for semi-diurnal terms 



,7^^ 



1 +fP multiplied by (l-2X)^ 1, {l + 2\f, for diurnal terms 
S 



.01' 



and l+^T^--^^ fo^^ the fortnightly term 

S 



(63) 



'* The following metliod of correcting the work of tlie paper on "Precession" has been rewritten, ancl 
was inserted on the 17th May, 1879. 
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Also we have the equations 



sill 4:€i 

sin 4e 

sin 2e\ 
sin4e 



1-X, 



sin 4e 
sin 4e 



1, 



sin 4€r 

^ 

sin 4e 



t /-. ^\ \ sin A€ 

i(i-2\), . 



sm 



_^ -^ sin 26^2 
4e ^' sin 46 



1 ~}~A 



i(i+2X). 



sin 4e'' 
sin 4e 



>^. . , (64) 



=\ 



Now we shall obtain a sufficiently accurate result, if the corrections be only applied 
to those terms in the differential equations which do not involve powers of q (or sin ^i)^ 
higher than the first. Then for the purpose of correction the differential equations to 
be corrected are by {77), (78), and (79) of Section 17 of ^^ Precession/' viz. :— ■ 



di 



ni' 



clt N ^n^ 



\^^p^q sin 4€^ + 2P^?(P^+3g'^) sin 2€\---|pg'(p^--c/^)^ sin 2e~\ 



clt 



L — rp^ sin. 4e]^/x 



■n 



dt 



y 



(65) 



As we are treating the obliquity as small, we may put 

\p^q=^p\{lo^+?>q^) = lpq{p^-(ff=lPQ and p^=P, 

when Pr=::cos i, Q=: sin i. 

Then for the purpose of correction, the terms depending on the moon's influence are 



cli 



ni' 



1 



clt N §7^, 



-^P^{sin 4€i+ sin 2^\ — sin ^e} 



dM 



> 



ni' 



clt 



2~-^sin4ei=::/x^^^ 



Uc 



• » • 



(66) 



And by symmetry- (or by (81) " Precession ") we have for the solar terms 



dt N^n* ^ 



dN',r- 



ni/ ■ 



clt 



:± ' p sin 4€ , 



« « 



• • {^'^) 



For the terms depending on the joint action of the sun and moon we have, by (82) 
and (33) '^ Precession,'' when the obliquity is treated as small, 



di 



'mm J 



TT 



dt 

clK 



Nq}i 



'iPQ sin 2e' 







> 



• • 



« * 



(68) 



f/im./ 



dt 
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Then if we multiply each of the sines by its appropriate factor given in (63), and 
substitute from (64) for each of them in terms of sin 4€, and collect the results from 
(66), (67), and (68), and express by the symbol § the corrections to be introduced for 
the effects of inertia^ we have 



^di 



/^ 



clt 
It 



^0 a 



N g^c * 



f_^—r 



h~~P .%^%%l-\f+Tf 



clt 



8^0 



n 



2 



Now 4{l-X)H5(l-2X)^— -| — (1-2X)(4-7X+4X2). Therefore if we add these 

d'i/ dlSJ" dp 

corrections to the full expressions for — , — (in which I put 1— |^^^=P) and /x—, given 



df dt 



dt 



w 



in (83) "Precession,'^ and write K=~7-5-~ for brevity, we have 



di 1 sin 4e 1 p^ j^2(i _ ^^ +^;_^^^4.K[(1 _2X) (1 _|X+X2)^+^;_l^^^)] 



dt N g% 



P 



^^l^Hi£pj^3(l_X)+T/ + |TT,^f + K[(l~X)V + T/] 



dt 



% 



S^^i 







p 



.3 



a--=^ — sm 4er 

' dt '^^nQ 



X 



l-^+K(l-^X)^ 



P 



^ 



(69) 



The last of these equations may be written approximately 



dt 
fjbd^ 



^ — sin 4€ P( 1 — ~ 



2 



8'^^. 







P; 



l^K(l-X)^]. 



• ' (70) 



Then if we multiply the two former of equations (69) by (70), and notice that, when 
P is taken as unity, 



2X' 



(1™2X)(1-^|X+X^)»- 1— p (l-^X)^:==iX(l--2X), 



and that 



1-(1-X)^=X(2~X) and ^i+(l-X)^=i(l-2X)(3-2X). 



we have 
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log tan^ 



ixd^ 



2 



1- 



2X /T, 



-' +K 



iX(l-2X) + X(2-X) -' +Jr(l-2X)(3 



2X) 



n(i 



clN 

1^4 



~^ + C-)\^^(}) + KX(2-X)\ 



1- 



X 
P 



T/J 



X 

T 



) 



J 



(71) 



If K be put equal to zero, we have tlie equations (84) which, were the subject of 
integration in Section 17 '^Precession/' 

Since K, \ and rf-^r^ are all small, the correction to the second equation is 
obviously insignificant, and we may take the term in K in the numerator of the first 
equation as being equal to ^K(l — 2X)(3 — 2X)(t^-^t). This correction is small although 
not insensible. This shows that the amount of change of obliquity has been slightly 
under-estimated. It does not, however, seem worth while to compute the corrected 
value for the change of obliquity in the integrations of the preceding paper. 

The equation of conservation of moment of momentum, which is derived from the 
integration of the second of (71), clearly remains sensibly unaffected. 

We see also from (70) that the time required for the changes has been over- 
estimated. If Kq, Xq ; K, X be the initial and final values of K and X at the beginning 
and end of one of the periods of integration ; then it is obvious that our estimate of 
time should have been multiplied by some fraction lying between 1 — Kq(1"— Xq)^ and 
1_K(1""X)^. 

Now at the beginning of the first period Kq='0364 and Xq==*0365, and at the end 
K=-0865 and X:=-0346. 

Whence Ko(l -Xo)^== '034, K(l -X)^=: -080. 

Hence it follows that the time, in the first period of the integration of Section 15, 
may have been over-estimated by some percentage less than some number lying 
between 3 and 8. 

In fact, I have corrected the first period of that integration by a rather more tedious 
process than that here exhibited, and I found that the time was over-estimated by a 
little less than 3 per cent. And it was found that we ought to subtract from the 
46,300,000 years comprised within the first period about 1,300,000 years. I also found 
that the error in the final value of the obliquity could hardly amount to more than V or 2". 

In the later periods of integration the error in the time would no doubt be a little 
larger fraction of the time comprised within each period, but as it is not interesting to 
find the time in anything but round numbers, it is not worth while to find the 
corrections. 

There is another point worth noticing. It might be suspected that when we 
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approach the critical point where n co&i=2I2, where the rate of change of obliquity 
was found to vanish, the tidal movements might have become so rapid as seriously to 
affect the correctness of the tidal theory used ; and accordingly it might be thought 
that the critical point was not reached even approximately when n cos i=2Ii, 

The preceding analysis will show at once that this is not the case. Near the critical 
point the solar terms ha^e become negligeable ; then if we put t^=0 in the first of 
equations (69) we have 

f =i-^/^" 4c.iP^[l-2\seci+K(l-2X)(l-|X+X^)]. 

The condition for the critical point in the first approximation was 2Xsec$=l ; if 
then i is so small that we may take sec^=l in the inertia term, this condition also 
causes the inertia term to vanish. 

Hence the corrected theory of tides makes no sensible difference in the critical point 

where -7: changes sis^n. 

Having now disposed of these special points connected with previous results, I 
shall return to questions of general dynamics connected with the approximate solution 
of the forced vibrations of viscous spheroids ; that is to say, I shall compare the results 
with those of — 

The forced oscillations of Jiiiid spheroids,'^ 

The same notation as before will serve again, and the equations of motion are 



dp . dW da 



UjJU (XtAj Cvy 

two similar equations 

and. —+----)----=: 

dx dy dz 



(73) 



If the external tide-generating forces be those due to a potential per unit volume 
equal to ivr^'^i, and r=a+cr^ be the equation to the tidal spheroid, where Si, di are 
surface harmonics of the i^^ order, then we must put 



W=t.(r^S,+2^"f -JM a-,+(3a^-r^)^ 



the second term being the potential of the tidal protuberance, and the last of the mean 
sphere. 

Differentiate the three equations of motion by x, y, % and add them, and we have 



V^(p"— ?^(3a^— 7 



OA 9 



^>V--o- 



^ This is a slight modification of Sir W. Thomsoh's investigation of the free oscillations of fluid 
spheres, Phil. Traos., 1863, p. 608. 
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j9=«^(3a^— r^)r~+ solid harmonics + a constant. 



Now when r=:a, at the mean surface of the sphere, p-i^gwo-i, therefore 



p=^(c^a--.^3)^^+^^^y . 



Then substituting this value of j^ in the equations of motion (73^, 



da 


d 


dt' 


dx 


da 


d 


dt' 


dx 



^^i+iir&+^^-'^t-^^^^^ 



2a 



OA fJ 



•g^A 



r 



2a ^^\a, 



' .^ 2(i-l) /rV 
and two similar equations 



(74) 



The expression within brackets [ ] on the right is the effective disturbing potential, 
inclusive of the effects of mutual gravitation, and thus this process is exactly parallel 
to that adopted above in order to include the effects of mutual gravitation in the dis- 
turbing potential in the case of the viscous spheroid. 

Now p, the radial velocity of flow, is equal to a-+/S-+y^, 

Therefore multiplying the equations (74) by -, -, - and adding them, we have, by 
the properties of homogeneous functions, 



But when r=a, p= 
Therefore 





dp 
dt 


ddi 




dt' 






d^CTi 



4-1 



S,:— 



2(^-1) 7^ ' 
2i + l9 d^' 



df 



iw 



■iq _ 2%{i-l) g 
^' 2i-\-l a 



(Ti, 



» > 



. . . (75) 



Now suppose Si=Q.i cos vt, and that the tidal motion is steady, so that di must be of 
the form XQi cos vt ; then substituting in (75) this form of cr/, we find 






'—V'^ + 



, . ^i(i-^ 9 



2i+l a 



•la 



i-i 



Whence 



id"^ 



(ri= 



2i(i^l) g 



■Qi cos vl 



(76) 



2^ + 1 a 



-—■'V' 



This gives the equation to the tidal spheroid. 
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Since the equilibrium tide, due to the disturbing potential, would be given by 



a' 



CT; — 



2 (^- 1) g 

2^ + l a 



Qi COS Vt, 



it follows that inertia augments the height of tide in the proportion 1 : 1—-^. — :~. —v^. 



V 



In the case where i=^2, the augmentation is in the proportion 1 : 1 — |— . 

8 

We will now consider the nature of the motion by which each particle assumes its 
successive positions. 

With the value of cri given in (76) 



2^ 4- 1 c(} 



— Qi'?^' 



2i{i~l) g 



cos vt. 



(2^ + l) a 



^■—v 



Then substituting in (74) 



da 



d v^ cos vt Qir^ 
dx 2i{i — \) g 



■v 



2i-\-l a 
and two similar equations 



V 



J 



Integrating with regard to t 



(77) 



d Ojrh sin vt 






^m 



^^ ^(tzD i'^^s 



2^+l a 
and two similar equations 



r 



(78) 



There might be a term introduced by integration, independent of the time, but this 
term must be zero, because if there were no disturbing force there w^ould be no flow. 
Hence it is clear that there is a velocity potential 9-, and that 



:^= 



d 



2i-\-l a 



(r^S.) 



(79) 



■V 



Now however slowly the motion takes place, there will always be a velocity potential, 
and if it be slow enough we may omit v^ in the denominator of (79). In other words, 
if inertia be neglected the velocity potential is 
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For the sake of comparison with the approximate solution for the tides of a viscous 
spheroid^ a precisely parallel process will now be carried out with regard to the fluid 
spheroid. 

We obtain a first approximation for -r ^ when inertia is neglected, by omitting v^ in 

ctv 

the denominator of (77); whence 



dot. 



cl I 2^4-1 g 

J — ^ 'I. 

clx\2i{i—l) a 



v^ ooB vtr'Qi). 



Substituting this approximate value in the equations of motion (73) we have 



dp , - d f^^y , 2i + l a r, 
J +7 1 W+t(^^T7r-"7. -t^-^ cos vtr .^, 
cldo dx\ 2^(^ — 1) a 

and two similar equations 







> 



(80) 



From these equations it is obvious that the second approximation to the form of the 

tidal spheroid is found by augmenting the equilibrium tide due to the tide-generating 

2^+1 a 
potential r^Q^ cos vt in the proportion 1 + Tpr. — 7x ""^^ to unity^ 

o 

When ^=2 the augmenting factor is 1+-|— . 



This is of course only an approximate result;. the accurate value of the factor is 



v 






1^ 1-^1— \ and we see that the two agree if the squares and higher powers of | 

are neghgeable. 

Now in the case of the viscous tides we found the augmenting factor to be 

1+~A%~' cos^ ^- When €::==0j which corresponds to the case of fluidityj the expressions 



are closely alike, but we should expect that the 79 ought really to be 75. 

The explanation which lies at the bottom of this curious discrepancy will be most 
easily obtained by considering the special case of a lunar semi-diurnal tide. 

We found in Part II., equation (21), the following values for a, /3, y. 



IDT 



^ 



""^Wv ^'" ^^ [{8a^--5r'');(/+4a;2y] 



/3 = ^ sin 2 € [( 8 a^ — 5 r~) x + ixy'^) > . 



lor 



y— ~ sm 2eAxy: 



V 



» 6 » e 



• (81) 
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where 

x=^r sin d cos lfh--o)t) 

yzzzr sin sin (c^— coi^) }>. 
ZT-r cos ^ 

Now consider tlie cavSa when the viscosity is infinitely small: here e is small, and 



sin 2€=:tan 2€=: 



""Svoy 



o 



5Qwa- 



r 



COT » COT 

Hence — sin 2€=— — , which is independent of the viscosity. 

By substituting this value in (81), we see that however small the viscosity, the 
nature of the motion, by which each particle assumes its successive positions, always 
preserves the same character ; and the motion always involves molecular rotation. 

But it has been already proved that, however slow the tidal motion of a fluid 
spheroid may be, yet the fluid motion is always irrotational. 

Hence in the two methods of attacking the same problem, different first approxi- 
mations have been used, whence follows the discrepancy of 79 instead of 75. 

The fact is that in using the equations of flow of a viscous fluid, and neglecting 

inertia to obtain a first approximation, we postulate that iv-z-, %o-j- tv-j-, are less im- 

•^•*- 5 r dt clt clt 

portant than t; V ^a, v V % vV^y\ and this is no longer the case if v be very small. 

It does not follow therefore that, in approaching the problem of fluidity from the 
side of viscosity, we must necessarily obtain even an approximate result. 

But the comparison which has just been made, shows that as regards the form of 
the tidal spheroid the two methods lead to closely similar results. 

It follows therefore that, in questions regarding merely the form of the spheroid, 
and not the mode of internal motion, we only incur a very small error by using the 
limiting case when t;==0 to give the solution for pure fluidity. 

In the paper on '^ Precession '' (Section 7), some doubt was expressed as to the 
applicability of the analysis, which gave the effects of tides on the precession of a 
rotating spheroid, to the limiting case of fluidity ; but the present results seem to 
justify the conclusions there drawn. 

The next point to be considered is the effects of inertia in — 

The forced oscillations of an elastic sphere. 

Sir William Thomson has found the form into which a homogeneous elastic sphere 
becomes distorted Uiuder the influence of a potential expressible as a solid harmonic 
of the points within the sphere. He afterwards supposed the sphere to possess the 
power of gravitation, and considered the effects by a synthetical method. The result 
is the equilibrium theory of the tides of an elastic sphere. When, however, the 
disturbing potential is periodic in time this theory is no longer accurate. 

4 F 2 
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It has already been remarked that the approximate solution of the problem of 
determining the state of internal flow of a viscous spheroid when inertia is neglected, 
is identical in form with that which gives the state of internal strain of an elastic 
sphere ; the velocities a, 13, y have merely to be read as displacements^ and the 
coefficient of viscosity v as that of rigidity. 

The effects of mutual gravitation may also be introduced in. both problems by the 
same artifice ; for in both cases we may take^ instead of the external disturbing 

potential wi^'^^Soosvi. an effective potential wr^{ScoBvt'—^~l and then cieem the 

sphere free of gravitational power. 

Now Sir William Thomson's solution shows that the surface radial displacement 
(which is of course equal to cr) is equal to 

Yq--( b cost^^— g-1 . . (82) 

1 Q Q 

If therefore we put (with Sir William Thomson) r^-^, we have -= cos^y^ 

^ ^ ' Zwcf? a r + g 

This expression gives the equilibrium elastic tide, the suffix being added to the 

cr to indicate that it is only a first approximation. 

Before going further we may remark that 

Scost;*^— ■tt~"="~— -Scos'?;^ ..,,,,,, (83) 

When, we wish to proceed to a second approximation^ including the effects of 
inertia, it must be noticed that the equations of motion in the two problems only 
differ in the fact that in that relating to viscosity the terms introduced by inertia are 

da d0 ^7 1 -1 , • ii r^ 1 • •, ,T (^^^'^ d^0 d^y 

•—lo-T'. -—w-^. _— whilst m the case oi elasticitv they are --iv—z, —iv-r^, —w~:r^. 

dt' dt' dt' J J clt^^ ^^3' ^ip. 

Hence a very slight alteration will make the whole of the above investigation 
applicable to the case of elasticity ; we have, in fact, merely to differentiate the 
approximate values for a, ^, y twice with regard to the time instead of once. 

Then just as before, we find the surface radial displacement, as fa,r as it is due to 

inertia, to be (compare (55)) 

wv^a? 79 Y 

cosi^^, 



9 o 



v^ 2.3.19^^7 
Y . . a- 

and— g cos t;^ must be put equal to (the first approximation) Scos'?;^— g-. Hence by 

tuhhi^ 79 X 
(57) and (83) the surface radial displacement due to inertia is — ^ 9 i q2 ~ — ^ ^^^ '^^' 

To this we must add the displacement due directly to the effective disturbing 

potential wr'^fS cos vt — ^-), where cr is now the second approximation, This we 

know from (82) is equal to 
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r— lb COS 1?^ — g - 1. 



Hence the total radial displacement is 

-TT— I b COS Vt'—U-+~77r"'^T?7: ' ^ COS vt 

19u \ ^a 19u 150 r + S 

But the total radial displacement is itself equal to cr. 
Therefore 

t- rr: S COS vt — U" + 7VX7 ^S COS Vt, 

a "^a ' 150(r-f g) 



nd. 

\ "^150rr + cry* 



o- S /, , 79-2;3 

- = -~7-- cos 'yn 1 +7^777 — X 
a t + g \ 150 (r + g} 



This is the second approximation to the form of the tidal spheroid, and from it we 
see that inertia has the effect of increasing the ellipticity of the spheroid in the 

proportion I+tttt; :■. 

^ ^ 160(r + g) 

Analogy with (76) would lead one to believe that the period of the gravest vibration 

/ 79 \"^ 
of an elastic sphere is 2wij^]; this result might be tested experimentally. 

If Q be put equal to zero, the sphere is devoid of gravitation, and if X be put 
equal to zero the sphere becomes perfectly fluid ; but the solution is then open to 
objections similar to those considered, when viscosity graduates into fluidity. 

It is obvious that the whole of this present part might 1)e easily adapted to that 
hypothesis of elastico- viscosity which was considered in the paper on '' Tides," but it 
does not at present seem worth while to do so. 

By substituting these second approximations in the equations of motion again, we 
might proceed to a third approximation, and so on ; but the analytical labour of the 
process would become very great. 

IV. Discihssion of the applicability of the results to the history of the earth. 

The first paper of this series was devoted to the consideration of inequalities of short 
period, in the state of flow of the interior, and in the form of surface, produced in a 
rotating viscous sphere by the attraction of an external disturbing body : this was the 
theory of tides. The investigation was admitted to be approximate from two causes 
— (i) the neglect of the inertia of the relative motion of the parts of the spheroid ; 
(ii) the neglect of tangential action between the surface of the mean sphere and the 
tidal protuberances. 
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In the second paper the inertia was still neglected^ but the effects of these tan- 
gential actions were considered, in as far as they modified the rotation of the spheroid 
as a whole. In that paper the sphere was treated as though it v/ere rigid, but had 
rigidly attached to its surface certain inequalities, which varied in distribution from 
instant to instant according to the tidal theory. 

In order to justify this assumption, it is now necessary to examine whether the tidal 
protuberances may be regarded as instantaneously and rigidly connected with the 
rotating sphere. If there is a secular distortion of the spheroid in excess of the 
regular tidal flux and reflux, the assumption is not rigorously exact ; but if the dis- 
tortion be very slow, the departure from exactness may be regarded as insensible. 

The first problem in the present paper is the investigation of the amount of secular 
distortion, and it is treated only in the simple case of a single disturbing body, or 
moon, moving in the equator of the tidally-distorted spheroid or earth. 

It is found, then, that the form of the lagging tide in the earth is not such that the 
pull, exercised by the moon on it, can retard the earth's rotation exactly as though the 
earth were a rigid body. In other words, there is an unequal distribution of the tidal 
frictional couple in various latitudes. 

We may see in a general way that the tidal protuberance is principally equatorial, 
and that accordingly the moon tends to retard the diurnal rotation of the equatorial 
portions of the sphere more rapidly than that of the polar regions. Hence the polar 
regions tend to outstrip the equator, and there is a slow motion from west to east 
relatively to the equator. 

When, however, we come to examine numerically the amount of this screwing 
motion of the earth's mass, it appears that the distortion is exceedingly slow, and 
accordingly the assumption of the instantaneous rigid connexion of the tidal protube- 
rance with the mean sphere is sufficiently accurate to allow all the results of the paper 
on ^^ Precession '^ to hold good. 

In the special case, which was the subject of numerical solution in that paper, we 
were dealing with a viscous mass which in ordinary parlance would be called a solid, 
and it was maintained that the results might possibly be applicable to the earth within 
the limits of geological history. 

Now the present investigation shows that if we look back 45,000,000 years from 
the present state of things, we might find a point in lat. 30° further west with reference 
to a point on the equator, by 4f^ than at present, and a point in lat. 60° further west 
by 14^^ The amount of distortion of the surface strata is also shown to be exceedingly 
minute. 

From these results we may conclude that this cause has had little or nothing to do 
with the observed crumpling of strata, at least within recent geological times. 

If, however, the views maintained in the paper on *^ Precession ^' as to the remote 
history of the earth are correct, it would not follow, from what has been stated above, 
that this cause has never played an important part; for the rate of the screwing of the 
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eartli s mass varies inversely as the sixth power of the moon's distance, multiplied by 
the angular velocity of the earth relatively to the moon. And according to that 
theory, in very early times the moon was very near the earth, whilst the relative 
angular velocity was comparatively great. Hence the screwing action may have been 
once sensible.'''^ 

Now this sort of motion, acting on a mass which is not perfectly homogeneous, 
would raise wrinkles on the surface which would run in directions perpendicular to the 
axis of greatest pressure. 

In the case of the earth the wrinkles would run north and south at the equator, and 
would bear away to the eastward in northerly and southerly latitudes ; so that at the 
north pole the trend would be north-east, and at the south pole north-west. Also the 
intensity of the wrinkling force varies as the square of the cosine of the latitude, and 
is thus greatest at the equator, and zero at the poles. Any wrinkle when once formed 
would have a tendency to turn slightly, so as to become more nearly east and west, 
than it was when first made. 

The general configuration of the continents (the large wrinkles), on the earth's sur- 
face appears to me remarkable when viewed in connexion with these results. 

There can be little doubt that, on the whole, the highest mountains are equatorial, 
and that the general trend of the great continents is north and south in those regions. 
The theoretical directions of coast line are not so well marked in parts removed from 
the equator. 

* This result is not strictly applicable to the case of infinitely small viscosity, because it gives a finite 

tliougli very small circulation, if tbe coefficient of viscosity be put equal to zero. 

By putting e=0 in (17'), Part I., we find a superior limit to the rate of distortion. With the present' 

dJj 
ans:ular velocities of the earth and moon, -:,- must be less than 5 X 10~^ cos^ in degrees per annum. 

^ clt 

dJj 
It is easv to find when — would be a maximum in the course of development considered in '* Preces- 

•^ dt 

sion;" for, neglecting the solar effects, it will be greatest when t^(w—Q) is greatest. 

Fow T^(^^— O) varies as [l+/<^— /*f — ^.|^~^]f~^^ and this function is a maximum when 

Taking /t=4-0074, andi=27-32, we have r^^l09-45r^ + 80-293=0. 

The solution of this is f=*2218. 

With this solution -— will be found to bo 56 million times as great as at present, being equal to 

clt 

18' cos^ per annum. With this value of ^, the length of the day is 5 hours 50 minutes, and of the 
month 7 hours 10 minutes. 

This gives a superior limit to the greatest rate of distortion which can ever have occurred. 

By (19'), however, we see that the rate of distortion per unit increment of the moon's distance may be 
made as large as we please by taking the coefficient of viscosity small enough. 

These considerations seem to show that there is no reason why this screwing action of the earth should 
not once have had considerable effects. (Added October 15, 1879.) 
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The great line of coast running from North Africa by Spain to Norway has a 
decidedly north-easterly bearing, and the long Chinese coast exhibits a similar ten- 
dency. The same may be observed in the line from Greenland down to the Gulf of 
Mexico, but here we meet with a very unfavourable case in Panama^ Mexico, and the 
long Californian coast line. 

From the paucity of land in the southern hemisphere the indications are not bo 
good, nor are they very favourable to these views. The great line of elevation which 
runs from Borneo through Queensland to New Zealand might perhaps be taken as an 
example of north-westerly trend. The Cordilleras run very nearly north and south, 
but exhibit a clear north-westerly twist in Tierra del Fuego, and there is another 
slight bend of the same character in Bolivia. 

But if this cause w^as that which principally determined the direction of terrestrial 
inequalities, then the view must be held that the general position of the continents 
has always been somewhat as at present, and that, after the wrinkles were formed, the 
surface attained a considerable rigidity, so that the inequalities could not entirely 
subside during the continu.ous adjustment to the form of equilibrium of the earth, 
adapted at each period to the lengthening day. With respect to this point, it is 
worthy of remark that many geologists are of opinion that the great continents have 
always been more or less in their present positions. 

An inspection of Professor Schiapparelli's map of Mars,'" I think, will prove that 
the north and south trend of continents is not something peculiar to the earth. In 
the equatorial regions we there observe a great many very large islands, separated by 
about tw^enty narrow channels running approximately north and south. The northern 
hemisphere is not given beyond lat. 40°, but the coast lines of the southern hemisphere 
exhibit a strongly marked north-westerly tendency. It must be confessed, however, 
that the case of Mars is almost too favourable, because we have to suppose, according 
to the theory, that its distortion is due to the sun, from which the planet must always 
have been distant. The very short period of the inner satellite shows, however, that 
the Martian rotation must have been (according to the theory) largely retarded ; and 
where there has been retardation, there must have been internal distortion. 

The second problem which is considered in the first part of the present paper is 
concerned with certain secondary tides. My attention was called to these tides by 
some remarks of Dr. Jules CAiiiiET,i" who says : — - 

" Les actions perturbatrices du soleil et de la lune, qui produisent les mouvements 
coniques de la precession des equinoxes et de la nutation, n'agissent que sur cette 
portion de Tellipsoide terrestre qui excede la sphere tangente aux deux poles, 
c'est-a-dire, en admettant Tetat p4teux de Tinterieur, a peu pres uniquement sur ce 

* ' Appendice alle Memorie deila Societa degli Spettroscopisti Italian!,' 1878, vol. vii., for a copy of 
wliich. I liave to thank M. Schiapparelli. 

t Societe Savoisienne d'Histoire et d'Archeologie, May 23, 1878. He is also autlior of a work, * Le 
Beplacement Pol aire.' I tliink Dr. Garret lias misunderstood Mr. Eyans. 
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que Ton est convenu d'appeler la croute terrestre^ et presque sur toute la croute tei^- 
restre. La croute glisse sur rinterieur plastique. EUe parvient k entrainer rinterieur, 
car, sinon, Taxe de la rotation du globe demeurerait parallele a Ini-meme dans Tespace, 
ou n eprouverait que des variations insignifiantes, et le phenomene de la precession des 
equinoxes n^existerait pas. Ainsi la croAte et Tinterieur se meuvent de quantites 
inegales^ d/ou le deplacement geographiqne du pole sur la sphere. 

^^ Cette idee a ete emise, je crois^ pour la premiere fois, par M. Evans ; depuis par 
M. J. Pekoche." 

Now with respect to this view, it appears to me to be sufficient to remark that, as 
the axes of the precessional and nntational couples are fixed relatively to the moon, 
whilst the earth rotates, therefore the tendency of any particular part of the crust to 
sHde over the interior is reversed in direction every twelve lunar hours, and therefore 
the result is not a secular displacement of the crust, but a small tidal distortion. 

As, however, it was just possible that this general method of regarding the subject 
overlooked some residual tendency to secular distortion, I have given the subject a 
more careful consideration. From this it appears that there is no other tendency to 
distortion besides that arising out of tidal friction, which has just been discussed. It 
is also found that the secondary tides must be very small compared with the primary 
ones ; with the present angular velocity of diurnal rotation, probably not so much in 
height as one-hundredth of the primary lunar semi-diurnal bodily tide. 

It seems out of the question that any heterogeneity of viscosity could alter this 
result, and therefore it may, I think, be safely asserted that any sliding of the crust 
over the interior is impossible — at least as arising from this set of causes. 

The second part of the paper is an investigation of the amount of work done in the 
interior of the viscous sphere by the bodily tidal distortion. 

According to the principles of energy, the work done on any element makes itself 
manifest in the form of heat. The whole work which is done on the system in a given 
time is equal to the whole energy lost to the system in the same time. From this 
consideration an estimate was given, in the paper on ^* Precession,^^ of the whole amount 
of heat generated in the earth in a given time. In the present paper the case is 
taken of a moon moving round the earth in the plane of the equator, and the work 
done on each element of the interior is found. The work done on the whole earth is 
found by summing up the work on each element, and it appears that the work per 
unit time is equal to the tidal frictional couple multiplied by the relative angular 
velocity of the two bodies. This remarkably simple law results from a complex law of 
internal distribution of work, and its identity with the law found in '* Precession," 
from simple considerations of energy, affords a valuable confirmation of the complete 
consistency of the theory of tides with itself 

Fig. 2 gives a graphical illustration of the distribution in the interior of the work 
done, or of the heat generated, which amounts to the same thing. The reader is 
referred to Part II. for an explanation of the figure. Mere inspection of the figure 

MDCCCLXXIX. 4 a 
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shows that by far tlie larger part of the heat is generated in the central parts, and 
calculation shows that abont one-third of the whole heat is generated within the central 
one-eighth of the volume, whilst in a spheroid of the size of the earth only one-tenth 
is generated within 500 miles of the surface. 

In the paper on "Precession^' the changes in the system of the sun, moon, and 
earth were traced backwards from the present lengths of day and month back to a 
common length of day and month of 5 hours 36 minutes, and it was found that in 
such a change heat enough must have been generated within the earth to raise its 
whole mass 3000^ Fahr. if applied all at once, supposing the earth to have the specific 
heat of iron. It appeared to me at that time that, unless these changes took place at 
a time very long antecedent to geological history, then this enormous amount of in- 
ternal heat generated would serve in part to explain the increase of temperature in 
mines and borings. Sir William Thomson, however, pointed out to me that the 
distribution of heat-o:eneration would probably be such as to prevent the realisation of 
my expectations. I accordingly made the further calculations, connected with the 
secular cooling of the earth, comprised in the latter portion of Part II. 

It is first shown that, taking certain average values for the increase of underground 
temperature and for the conductivity of the earth, then the earth (considered homo- 
geneous) must be losing by conduction outwards an amount of energy equal to its 
present kinetic energy of rotation in about 262 million years. 

It is next shown that in the passage of the system from a day of 5- hours 40 minutes 
to one of 24 hours, there is lost to the system an amount of energy equal to 13|- times 
the present kinetic energy of rotation of the earth. Thus it appears that, at the 
present rate of loss^ the internal friction gives a supply of hea,t for 3,560 million years. 
So far it would seem that internal friction might be a powerful factor in the secular 
cooling of the earth, and the next investigation is directly concerned with that question. 

In the case of the tidally-distorted sphere the distribution of heat -generation 
depends on latitude as well as depth from the surface, but the average law of heat- 
generation, as dependent on depth alone, may easily be found. Suppose, then, that 
we imagine an infinite slab of rock 8,000 miles thick, and that we liken the medial 
plane to the earth's centre and suppose the heat to be generated uniformly in time, 
according to the average law above referred to. Then conceive the two faces of the 
slab to be always kept at the same constant temperature^ and that initially, when 
the heat-generation begins, the whole slab is at this same tempera.ture. The problem 
then is, to find the rate of increase of temperature going inwards from either face of 
the slab after any time. 

This problem is solved, and by certain considerations (for which the reader is referred 
back) is made to give results which must agree pretty closely with the temperature 
gradient at the surface of an earth in which 13i times the present kinetic energy of 
earth's rotation, estimated as heat, is uniformly generated in time, with the average 
space distribution referred to. It appears that at the end of the heat-generation the 
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temperature gradient at the surface is sensibly the same, at whatever rate the heat is 
generated, provided it is all generated within 1,000 million years ; but the temperature 
gradient can never be quite so steep as if the whole heat were generated instan- 
taneously. The gradient, if the changes take place within 1,000 million years, is 
found to be about 1° Fahr. in 2,600 feet. Now the actually observed' increase of 
underground temperature is something like 1° Fahr. in 50 feet; it therefore appears 
that perhaps one-fiftieth of the present increase of underground temperature may pos- 
sibly be referred to the effects of long past internal friction. It follows, therefore, that 
Sir William Thomson's investigation of the secular cooling of the earth is not 
sensibly affected by these considerations. 

If at any time in the future we should attain to an accurate knowledge of the 
increase of underground temperature, it is just within the bounds of possibility that a 
smaller rate of increase of temperature may be observed in the equatorial regions than 
elsewhere, because the curve of equal heat generation, which at the equator is nearly 
500 miles below the surface, actually reaches the surface at the pole. 

The last problem here treated is concerned with the effects of inertia on the tides of 
a viscous spheroid. As this part will be only valuable to those who are interested in 
the actual theory of tides, it may here be dismissed in a few words. The theory used 
in the two former papers, and in the first two parts of the present one, was founded 
on the neglect of inertia ; and although it was shown in the paper on '^ Tides'' that 
the error in the results could not be important, in the case of a sphere disturbed by 
tides of a frequency equal to the present hmar and solar tides, yet this neglect left a 
defect in the theory which it was desu^able to supply. Moreover it was possible that, 
when the frequency of the tides was much more rapid than at present (as was found 
to have been the case in the paper on ''Precession"), the theory used might be 
seriously at fault. 

It is here shown (see (62)) that for a given lag of tide the height of tide is a little 
greater, and that lor a given frequency of tide the lag is a little greater than the 
approximate theory supposed. 

A rough correction is then appHed to the numerical results given in the paper on 
'' Precession" for the secular changes in the configuration of the system ; it appears 
that the time occupied by the changes in the first solution (Section 15) is overstated 
by about one-fortieth part, but that all the other results, both in this solution and 
the other, are left practically unaffected. To the general reader, therefore, the value 
of this part of the paper simply lies in its confirmation of previous work. 

From a mathematical point of view, a comparison of the methods employed with 
those for finding the forced oscillations of fluid spheres is instructive. 

Lastly, the analytical investigation of the effects of inertia on the forced oscillations 
of a viscous sphere is found to be applicable, almost verbatim, to the same problem 
concerning an elastic sphere. The results are complementary to those of Sir William 
Thomson's statical theory of the tides of an elastic sphere. 
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